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modest  memory  size.  Examples  are  included  to  show  the  versatility  of  this 
approach  to  the  optimization  problem,  as  well  as  its  limitations. 
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Array  Optimization  Using  Subarrays 

THE  CONCEPT 

INTRODUCTION 

Optimization  of  the  element  excitations  of  discrete  antenna  arrays  is  a 
matter  of  definition  for  three  reasons.  First,  the  definition  of  optimality 
will  dictate  the  appropriate  mathematical  approach.  Seemingly  subtle 
changes  in  the  definition  of  optimality  can  radically  alter  the  applicable 
mathematical  methods.  Second,  element  excitations  that  are  optimal  in  one 
sense  are  unlikely  to  be  optimal  in  another  sense.  Two  sets  of  excitations, 
each  set  optimal  in  its  own  sense,  can  be  completely  different.  Third,  the 
definition  of  optimality  must  reflect  directly  upon  the  primary  design  goals 
for  the  array.  It  does  no  good  to  optimize  the  directivity  index  and  then 
complain  that  the  sidelobes  are  too  high  because  the  design  goal  of  low 
sidelobes  and  the  definition  of  optimality  (maximum  DI)  are  not  directly 
related. 

This  report  defines  and  uses  exclusively  the  concept  of  directivity 
index  with  beamwidth  control  (DIBC).  Several  advantages  as  well  as 
difficulties  inherent  in  this  definition  are  discussed.  The  primary 
difficulty  in  this  definition  is  the  requirement  of  large  computer  memories 
for  large  arrays.  A technique  employing  subarrays  of  the  full  array  in  a 
systematic  manner  is  shown  to  overcome  this  problem. 

The  optimization  procedure  in  this  report  is  applicable  when  the 
following  premises  obtain: 

1.  The  wavelength  X of  the  design  frequency  is  given  and  fixed. 

2.  The  number  n of  elements  in  the  array  is  fixed  and  all  the  element 
positions  (x|^,  yj^,  z^^),  k = l,...,n  are  known  and  fixed. 

3.  Individual  element  field  patterns  at  the  design  frequency  are 
completely  known. 

4.  The  ambient  noise  field  at  the  design  frequency  is  completely  known. 

5.  Element  interactions  can  be  ignored. 

6.  Element  excitations  can  be  phased  (i.e.,  complex). 


FIELD  PATTERNS  AND  COORDINATE  SYSTEM 

The  spherical  coordinate  system  of  figure  1 is  used  throughout  this 
report;  however,  a particular  direction  (9,4>)  will  be  specified  by  the 
direction  cosines 
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cos  ot  = sin  <j>  cos  9 
cos  B = sin  4i  sin  6 
cos  Y = cos  <|)  . 


(1) 


Figure  1.  The  Coordinate  System 


The  most  general  field  pattern  treated  here  is 


V(0,(}>) 


(2) 


where  R]^(0,<t))  is  the  phased  response  of  the  k-th  element,  and 

dk(9,<t>)  = cos  CO®  ^ ^ ^k  ^ ‘ 

Because  of  assumptions  1-6,  the  field  pattern  V(0,<t))  depends  solely  upon  the 
phased  (complex)  excitations  ai,...,an.  The  ambient  noise  field  N(0,ij)) 
enters  only  in  the  definition  of  optimal  excitations. 


1 

i 
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DIRECTIVITY  INDEX  WITH  BEAMWIDTH  CONTROL  (DIBC) 

The  antenna  designer  is  required  to  divide  the  set  of  all  directions, 
denoted  into  three  disjoint  regions: 

■wi  = Mainlobe  region 
^ = Sidelobc  region 
cl  = Ignored  region  = - )• 

This  division  of  directional  space  is  completely  arbitrary,  except  that 
neither*?1'i  nor  can  be  empty  sets,  whereas  cS  can  be  empty  if  desired. 

Once  a particular  choice  of  7>L,  ^ , and  cH  has  been  made,  the  following 
definition  of  optimality  is  used. 

Definition  1.  The  element  excitations  a^,...,a^  are  optimal  excita- 
tions for  a given  choice  of  regions  ^ , and  c|  , if  and  only  if  the 

ratio 


DIBC 


//inlN(e.(<))  V^(0.<|>)|  sin  4*  d4>d0 
//  |n(0,<j>)  v2(9,(J,)|  sin  <|>  d4)d9 


(4) 
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is  maximized.  Any  ratio  of  this  form  will  be  referred  to  as  a Directivity 
Index  with  Beamwidth  Control  (DIBC). 

Maximizing  Directivity  Index  (DI)  is  a limiting  case  of  maximizing  DIBC. 
To  see  this,  recall  that  for  a specified  direction  ^ maxi- 

mum if  the  ratio 


DI 


I ♦o>  *0)  I 

//jjl  N(9,  (}))  V^(6,  (f)  1 sin  4)  d(t)d0 


(5) 


is  maximized.  Now  let  the  ignored  region cX  be  empty,  let  the  mainlobe 
region“?»l  contain  (Qq*  letj?f  = ^2  • Then  excitations 

maximizing  DIBC  converge  to  excitations  which  maximize  DI  as  the  mainlobe 
regi on  shrinks  down  on  the  point  (Qq,  4>o)" 


We  have  defined  optimal  excitations  as  those  for  which  DIBC  is  maximized 
for  some  choice  of  regions*^  > j and  «Jl  . This  allows  a measure  of  control 
over  the  beamwidth  and  sidelobe  level.  By  varying  systematically  the  choice 
oiVl  and  , and  with  each  choice  maximizing  the  DIBC,  we  can  directly 
examine  the  tradeoff  between  beamwidth  and  sidelobe  level  for  the  particular 
array  at  hand.  The  engineer  can  then  select  those  excitations  which  best 
suit  his  needs.  Generally,  the  larger  the  mainlobe  region*!^  and  the  smaller 
the  sidelobe  regionj?^  (for  fixed  ignored  region^  ),  the  lower  the  overall 
sidelobe  level  and  the  greater  the  beamwidth.  However,  this  is  not  always 
the  case,  since  sidelobe  level  does  not  enter  directly  into  the  DIBC  ratio 
of  equation  (4).  Nothing  prevents  the  field  pattern  from  having  narrow  high 
amplitude  sidelobes,  since  such  sidelobes  would  contribute  little  to  the 
denominator  of  the  DIBC. 


Another  reason  for  maximizing  DIBC  is  simply  that  it  is  conceptually 
easy  to  do  so.  All  that  is  required  is  the  solution  of  an  eigenvalue/ 
eigenvector  problem  (see  Theorem  1),  and  problems  of  this  type  have  been 
extensively  studied  in  the  literature  (reference  1).  Numerically,  such 
problems  require  considerable  care.  Fortunately,  well-designed  computer 
programs  are  available  for  the  solution  of  eigenproblems  (references  2 and 
3).  With  the  use  of  these  routines,  the  solutions  of  the  eigenproblems 
encountered  in  the  antenna  problem  seem  to  be  numerically  stable.  This  is 
not  to  say  that  there  may  not  be  arrays  which  yield  numerically  unstable 
e i genproblems . 

A final  reason  for  maximizing  DIBC  is  more  esoteric.  In  the  process  of 
solving  the  required  eigenproblem,  all  the  eigenvalue/eigenvector  pairs  can 
be  computed,  not  merely  the  largest  one.  It  happens  that  the  field  patterns 
corresponding  to  the  lower  order  eigenvalues  have  some  interesting  features 
(see  example  1,  figures  5-8).  In  addition,  it  happens  that  the  larger 
eigenvalues  are  "close  together";  i.e.,  several  linearly  independent  sets  of 
excitations  exist  which  give  DIBC  values  that  lie  "close  together."  (For  an 
analogous  situation,  see  reference  4.)  What  this  means  in  the  antenna 
problem  is  that  without  sacrificing  antenna  performance  (as  measured  solely 
by  the  DIBC),  it  becomes  a simple  matter  to  examine  numerous  different  sets 
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of  excitations  with  the  aim  of  improving  some  completely  different  design 
goal  of  the  array.  This  will  not  be  discussed  further  in  this  report. 

It  must  be  mentioned  that  this  approach  to  the  array  optimization  problem 
does  not  attempt  to  address  several  issues  which  are  of  practical  interest. 
First,  this  approach  does  not  guarantee  that  the  array  performance  is  insen- 
sitive to  perturbations  in  the  optimum  excitations.  The  question  of  sensi- 
tivity to  excitation  perturbation  can  be  examined  only  after  the  optimum 
excitations  are  found.  Second,  this  approach  does  not  attempt  to  control 
the  efficiency  of  the  array.  In  other  words,  it  can  happen  that  the  optimal 
excitations  for  a particular  array  may  drive  certain  elements  at  their  maxi- 
mum allowed  levels  while  the  remaining  elements  are  hardly  driven  at  all,  so 
that  the  total  power  out  of  the  array  is  too  low  for  the  application.  This 
problem  is  common  to  all  amplitude  shaded  arrays  and  can  be  examined  after 
the  optimum  excitations  are  found.  Finally,  this  approach  to  array  optimiza- 
tion requires  that  the  optimal  excitations  be  allowed  to  be  complex  (i.e., 
phased),  and  this  can  require  certain  elements  to  absorb  power  from  the 
medium.  We  stress  that  this,  as  well  as  the  other  two  practical  matters, 
must  be  addressed  after  the  optimal  excitations  are  found. 


COMPUTER  STORAGE  PROBLEM 

The  primary  drawback  to  maximizing  DIBC  is  that  the  number  of  computer 
storage  locations  required  (using  the  program  in  appendix  B)  is  approximately 

Nt  = 6n2  + 16n  + 12,000  , (6) 

for  the  case  of  constant  ambient  noise  field  and  omnidirectional  elements, 
(table  1).  Since  the  total  requirement  will  grow  as  the  ambient  noise  field 
and/or  element  field  patterns  require  more  storage  to  compute,  it  appears 
that  the  direct  computation  of  optimal  excitations  for  any  array  of  100  or 
more  elements  requires  large,  main  frame  computers. 


Table  1.  Storage  Requirements  for  Direct 
Computation  of  Maximum  DIBC 


No.  of  Elements,  n 

Minimum  Storage  Required,  Nf 

15 

13,590 

30 

17,880 

60 

34,560 

100 

73,600 

150 

149,400 

200 

255,200 

The  storage  requirements  for  maximizing  DIBC  can  be  avoided.  A technique 
known  as  group  coordinate  relaxation  (reference  5)  gives  a method  which  can 
be  tailored  to  the  amount  of  computer  memory  available.  The  technique  is  an 
excellent  example  of  how  to  trade  off  computer  memory  for  computational 
speed.  The  more  memory  available,  the  faster  the  DIBC  can  be  maximized. 


t 

f 
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Group  coordinate  relaxation,  in  the  context  of  maximizing  DIBC,  is  simply 
stated.  Suppose  there  are  300  elements  in  the  array.  Make  any  initial  guess 
at  the  optimal  excitations.  Define  distinct  subarrays  of,  say,  50  elements 
each.  By  working  with  the  first  of  these  subarrays,  new  element  excitations 
are  computed  for  these  50  elements,  so  that  the  DIBC  of  the  entire  300  ele- 
ment array  is  increased.  Next,  new  excitations  are  computed  for  the  second 
subarray.  Cycling  through  all  six  subarrays  in  turn,  until  DIBC  for  the 
entire  300  element  array  cannot  be  increased  further  by  changing  the  excita- 
tions in  any  of  the  subarrays,  is  the  essence  of  group  coordinate  relaxation. 
The  method  can  be  proved  to  be  convergent.  The  rate  of  convergence  depends 
heavily  on  the  size  of  the  subarrays  used.  The  larger  the  subarrays,  the 
faster  the  convergence,  and  the  more  core  storage  required.  Thus,  core  stor- 
age is  traded  off  in  a direct  manner  for  the  convergence  rate,  and  hence  for 
computation  time.  In  addition,  each  step  of  the  group  coordinate  relaxation 
method  produces  new  excitations  which  increase  the  DIBC,  so  that  if  the  com- 
putations are  interrupted  for  any  reason: 

(1)  The  last  computed  excitations  are  better  than  any  of  the  excitations 
previously  computed, 

(2)  By  saving  the  last  computed  excitations,  the  computations  can  be 
resumed  without  any  significant  loss. 

If  ng  is  the  number  of  elements  in  a subarray  used  by  the  group  coordinate 
relaxation  process,  the  total  storage  required  (using  the  program  in  appendix 
B)  is  approximately 


N„  = 6n^  + 8(n  + n ) + 12,000  , (7)  ' 

R s s ’ 

for  the  case  of  constant  ambient  noise  field  and  omnidirectional  elements 
(table  2). 


Table  2.  Storage  Requirements  for  Maximizing  DIBC  of 
200-Element  Array  by  Group  Coordinate  Relaxation 


No.  of  Subarray  Elements,  ng 

Minimum  Storage  Required,  Nr 

15 

15,070 

30 

19,240 

60 

35,680 

100 

74,400 

150 

149,800 

200 

255,200 

ELABORATION  OF  THE  CONCEPT 


DIBC  AND  THE  EIGENPROBLEM 

T 

Let  the  vector  a = <a,,...,a  > be  the  vector  of  element  excitations  for 

1 n 

the  field  pattern  V(0,4))  given  by  equation  (2).  Then 
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JJ  |N(e,4>)  v2(e,(t.)|  sin  <t>  d4)d0 
= //  lN(e,4))l  lv(0,(»))l^  sin  d<t)d0 


= //  |N(e,(|))| 


n 

Z 

k=l 


a^R^(0,(|>)exp 


2iri 

X 


d,  (0 
k 


. Z a.R.(0,<}i)exp d.(0,(j))  isin  4)  d4)d0 
j = l J J J J\ 


n n 

= Z Z a.a.  J/  |n(0,(}))|  R,(0,(t))  R.(0,4)) 
k=l  j = l ^ J 


exp 


j^jdj(0,(|))  - d^(0,(t))  si 


in  4>  d(t)d0 


— T 

= a Ua, 


where  U is  an  n x n complex  matrix.  If  U = with  k denoting  the  row 

number  and  j denoting  the  column  number,  then 


u 


= //  |N(0,(t))i  Rj(0,4))R^(0,(l))exp|m[dj 


(0,<t>)  - d^(0,(j))]| 


sin  <j)  d<t>d0. 


(8) 


Clearly,  U is  a Hermitian  matrix  (i.e.,  U = U^),  since  it  is  obvious  that 
Ukj  ~ ^jk"  Also,  U is  positive  definite  since 

a'^Ua  = //  1n(0,4.)  V^(0,(())|  sin  4)  d4)d0  > 0 , f9) 

whenever  the  excitation  vector  a 0 (and  provided  the  mainlobe  region  ■h  is 
not  a set  of  measure  zero,  a pathological  coi;dition  that  is  not  encountered 
in  this  application).  Therefore,  for  every  mainlobe  region  the  matrix  U 
defined  in  equation  (8)  is  an  n x n positive  definite  Hermitian  matrix. 
Similarly, 


//  |n(0,4;)  V^(0,4>)|sin  4>  d4id0  = a'^Wa  , 
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where  W - is  an  n x n positive  definite  Hermitian  matrix  whose  gen- 

eral entry  is 


"kj  RjO,*!))  R^(9,4>)exp|m  jdj(e,(l))-d^^(e,<|))j|sin  (|>  d4id9.  (10) 


Thus,  for  a given  choice  ofT)^  » and«A.  , we  have 


DIBC 


-T 
a Ua 

-T 
a Wa 


(11) 


which  is  a ratio  of  positive  definite  Hermitian  forms.  Therefore,  optimal 
excitations  are  those  that  maximize  this  ratio  of  Hermitian  forms. 

The  mathematical  tools  for  handling  ratios  of  the  form  of  equation  (11) 
have  been  known  for  at  least  a century.  We  have  the  following  general  mathe- 
matical result. 


Theorem  1.  If  U and  W are  n x n Hermitian  matrices  and  W is  positive  defin- 
ite, then  the  eigenvalues  of  the  generalized  eigenproblem 

Uz  = yWz  (12) 

are  all  real.  Let  yi  > P2  ^ 
linearly  independent  vectors 

Uzk  = 


...  ^ yn  denote  these  eigenvalues.  Then 
Z]^,...,Zj|  can  be  found  which  satisfy 

>R“l»«»*>n  (13) 


and 


The  vectors  zj,...,Zn  are  called 
( 12 ) . Also,  we  have 


z^  Wz  if  k = j 

\ j \0  if  k j . 


the  eigenvectors  of  the  eigenproblem 


(14) 


max 

z)*0 


-T 
z Uz 


-T 
z Wz 


(15) 


and  this  maximum  is  attained  for  every  eigenvector  corresponding  to  yj,  and 


min 

z^O 


-T 
z Uz 

-T 
z Wz, 


(16) 
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and  this  minimum  is  attained  for  every  eigenvector  corresponding  to 
Finally,  if  1 £ k £ n,  then  for  any  constants  not  all  zero,  we 

have 


y^  i z^Uz  - \ ^ 


where  z = ajzj  + ...  + Oii^zi^.  - 

The  proofs  of  the  various  parts  of  this  theorem  can  be  found  in  numerous 
sources,  e.g.,  reference  1. 

For  the  immediate  purposes,  the  most  important  part  of  this  theorem  is 
equation  (15).  It  states  that  optimal  excitations  are  precisely  the  compo- 
nents of  any  eigenvector  corresponding  to  the  largest  eigenvalue  of  the 
generalized  eigenproblem  Uz  = yWz,  where  U and  W are  defined  by  equations 
(8)  and  (10). 

Theoretically,  Theorem  1 solves  the  problem  of  maximizing  DIBC  complete- 
ly. A discrete  reformulation  of  DIBC  is  discussed  next,  followed  by  a dis- 
cussion of  the  numerical  methods  for  the  solution  of  the  required  eigenprob- 
lem. 


A DISCRETE  VERSION  OF  DIBC 


Maximizing  the  DIBC  ratio  (4)  is  mathematically  tractable,  but  it  is  not 
practical.  It  requires  the  solution  of  an  eigenproblem,  which  in  turn 
requires  the  evaluation  of  approximately  n^  double  integrals  over  subsets 
of  the  unit  sphere.  Since  it  is  essential  that  the  mainlobe  region ^ and 
the  sidelobe  region  V* be  quite  general  in  nature  (i.e.,  be  defined  to  suit 
the  particular  situation),  these  double  integrals  are  impossible  to  evaluate 
explicitly  in  general,  and  are  also  difficult  and  time  consuming  to  evaluate 
accurately  by  numerical  methods.  For  these  reasons,  DIBC  itself  is  NOT  opti- 
mized. What  is  optimized  is  a discrete  version  (DIBCFp)  of  DIBC  that  is 
not  only  numerically  practical  to  use,  but  is  also  conceptually  simple. 


I 


The  discrete  DIBC  definition  replaces  the  surface  integrals  in  ratio  (4) 
by  discrete  sums  over  points  chosen  in"?^  an  Since'?*J_  and  ^ are  not  known 

a priori,  these  points  are  distributed  "uniformly"  over  the  surface  of  the 
sphere,  with  each  point  contributing  one  term  to  the  discrete  sum,  and  all 
terms  entering  with  equal  weight.  Ideally,  then,  these  points  must  show  no 
directional  bias  and  must  be  easy  to  compute.  Furthermore,  it  must  be  pos- 
sible to  choose  these  points  with  any  desired  "density"  on  the  sphere. 

A natural  choice  for  points  fulfilling  these  conditions  is  easy  to  de- 
scribe, but  difficult  to  compute:  Choose  as  points  the  equilibrium  posi- 
tions of  a finite  number  of  positive  charges  constrained  to  lie  on  the  sur- 
face of  the  unit  sphere.  When  the  number  of  positive  charges  is  4,  6,  8, 

12,  or  20,  it  is  intuitively  clear  that  stable  points  for  these  charges  are 
at  the  vertices  of  the  5 regular  Platonic  bodies:  the  tetrahedron,  the 
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octahedron,  the  cube,  the  icosahedron,  and  the  dodecahedron,  respectively. 
Unfortunately,  these  are  the  only  easy  cases  (see  reference  6). 

The  discrete  points  chosen  to  define  discrete  DISC  are  the  vertices  of  a 
geodesic  dome.  Consider  the  icosahedron  shown  in  figure  2. 


Figure  2.  The  Icosahedron 


(Note  that  in  this  figure  the  3r-axis  is  in  the  plane  of  the  paper  and  the 
z-axis  is  tilted  slightly  to  show  off  the  configuration.  The  x-axis  is  not 
shown,  but  is,  of  course,  orthogonal  to  the  yz-plane.)  This  regular  figure 
has  12  verti^s^  faces , ^d^^30,  ^ donjgs.  wLth  - • 

’num1?er*  ^ Taces  are  constructed  from  the  icosahedron  by  subdividing  its 
equilateral  triangular  faces  in  a systematic  manner  (reference  7).  First, 
subdivide  each  face  into  congruent  equilateral  subtriangles  as  shown  in 
figure  3;  that  is,  for  each  positive  integer  p ^ 1,  find  p+1  equispaced 
points  along  each  edge  and  pass  lines  through  each  of  these  points  parallel 
to  the  other  two  edges.  Next,  take  all  the  vertices  of  the  equilateral 
subtri angles  so  generated  and  project  them  onto  the  unit  sphere.  By  doing 
this  for  each  face  of  the  icosahedron  for  a fixed  integer  p > 1,  we 
construct  the  vertices  of  a geodesic  dome  of  order  p.  We  define  the  Fuller 
points  JJp  to  be  the  totality  of  these  points. 


Figure  3.  One  Face  of  Icosahedron  Subdivided  into  p Parts 
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I 

The  Fuller  points  are  uniquely  oriented  in  Cartesian  space  once  the 
vertices  of  the  icosahedron  are  defined.  With  some  simple  trigonometry, 
table  3 results.  How  many  points  are  there  in^p?  By  inspecting  an 


Table  3.  Vertices  of  the  Icosahedron 


X 


000000000 

894427191 

276393202 

723606798 

723606798 

276393202 

723606798 

276393202 

894427191 

276393202 

723606798 

000000000 


y 

. 000000000 
.000000000 
.850650808 
.$25731112 
-.525731112 
-.850650808 
.525731112 
.850650808 
.000000000 
-.850650808 
-.525731112 
.000000000 


1 . 000000000 
.447213595 
.447213595 
.447213595 
.447213595 
.447213595 

- .447213595 

- .447213595 

- .447213595 

- .447213595 

- .447213595 

-1.000000000 


"unfolded"  paper  model  of  the  icosahedron,  it  is  easy  to  see  that^p, 
contains  exactly  10  p^  + 2 points.  Table  4 gives  an  idea  of  the  density 
of  these  points  on  the  unit  sphere. 


T ebl ^4 « .Number  and  Density  of  Fuller  PointsJJp 


p 

No.  of  Points  in  ^p 

No.  of  Steradians/Point 

1 

12 

1.047 

2 

42 

.299 

4 

162 

.0776 

8 

642 

.0196 

16 

2562 

. 00490 

24 

5762 

.00218 

32 

10,242 

.00123 

Notice  that  the  Fuller  points ^p  are  not  quite  ideal.  Those  points 
chosen  near  the  center  of  a face  of  the  original  icosahedron  will  be  less 
finely  spaced  when  projected  onto  the  sphere  than  will  those  points  that 
were  chosen  nearer  an  edge.  This  defect  in  ^p  does  not  seem  to  be 
significant  in  this  application. 

With  the  Fuller  points  defined,  we  state  the  following. 

Definition  2.  For  a given  integer  p ^ 1,  and  regions7»J  , and 
the  element  excitations  a]^,...,a{,  are  optimal  if  and  only  if  the  ratio 
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DIBCF 

P 


(B  t^)c^ 

E 


1 E 

j(e  ,4>  )e  n (^  (^ ) 


|N(6,(j>)V^(e,(t))| 
|N(e,(j))  v^(e,(t))| 


(18) 


is  maximized,  where  is  the  number  of  points  of  ^ p that  lie  in*^ 

l«^pn(^U4)|  is  the  number  of  points  oftfp  that  lie  inn^j(  Any  ratio  of 
the  form  of  equation  (18)  will  be  referred  to  as  a Directivity  Index  with 
Beamwidth  Control  over  the  Fuller  points  (DIBCFp). 


The  formulation  of  DIBCFp  as  an  eigenproblem  parallels  that  for  DIBC. 
Specifically,  we  have 


Mz  = pSz  , (19) 

where  M = and  S = [sj^j^  are  n x n positive  definite  Herraitian 

matrices  with 


n,  . 

•tj 


n 'P(  T I (e,i(.)€j  n*^ 


N(e,()>)|Rj(e,(|i)R|^(9,(>)exp^2iti  [dj(e,4i)-d^{e,*)]^ 


(20) 


(6,<^)ca^nc\u/) 


I N(  6 , ♦ ) I R^  ( e , 4.  )R^TffT?ye  xp^m  [d  j ( e , )-d|^(  e , ♦ )]j 


(21) 


By  Theorem  1,  maximizing  DIBCFp  requires  the  computation  of  any  eigenvec- 
tor corresponding  to  the  largest  eigenvalue  for  the  eigenproblem  (19).  The 
numerical  solution  of  (19)  is  fully  discussed  in  the  next  section. 

There  are  two  considerations  that  should  enter  into  the  particular  choice 
of  p for  the  Fuller  points  *^p.  First,  the  points  inU'p  should  be  numer- 
ous enough  to  adequately  sample  the  worst  behavior  of  any  realizable  field 
pattern.  In  other  words,  p should  be  large  enough  that  even  the  narrowest 
sidelobe  achievable  in  the  field  pattern  will  contain  points  in  ^p. 

Second,  Theorem  1 requires  that  the  denominator  matrix  S of  the  DIBCFp 
ratio  be  positive  definite.  Normally,  the  sampling  criterion  will  effect 
this . 


NUMERICAL  SOLUTION  OF  THE  EIGENPROBLEM:  DIRECT  METHOD 

The  eigenproblem  (19)  is  equivalent  to  the  eigenproblem 

(S"l  M)z  “ yz  , (22) 

In  other  words,  the  eigenvalues  and  eigenvectors  of  equation  (22)  are  pre- 
cisely the  same  as  those  of  equation  (19).  There  are  two  difficulties  in 
using  equation  (22)  for  nxnerical  computation.  First,  it  requires  the 
inverse  of  the  matrix  S,  whose  only  special  structure  is  that  it  is  positive 
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detinite  and  Hermician.  In  general,  numerical  computation  of  the  inverse  of 
matrices  should  be  avoided  if  possible.  Second,  equation  (22)  is  not  a Her- 
mit i an  eigenprobl em;  that  is,  S~^M  is  not  necessarily  Hermitian  even  though 
S and  M are  both  Hermitian.  This  means  that  the  eigenvalues  and  eigenvectors 
of  equation  (22)  must  be  computed  by  a routine  designed  for  a general  complex 
matrix,  and  this  means  that  the  eigenvalues  can  (and  do)  turn  out  to  be  com- 
plex numbers  because  of  numerical  roundoff.  Since  Theorem  1 requires  that 
all  the  eigenvalues  be  strictly  real  numbers,  there  is  numerical  error  in 
using  equation  (22)  caused  by  destruction  of  the  natural  Hermitian  symmetry 
in  equation  (19).  For  these  reasons,  it  is  desirable  to  solve  the  eigenprob- 
I em  (19)  directly. 

Martin  and  Wilkinson  (reference  3)  contain  a method  and  a routine  for 
solving  this  eigenproblem  when  M and  S are  real  syranetric.  Both  the  tech- 
nique and  the  routine  can  be  adapted  to  the  Hermitian  case.  Every  Hermitian 
positive  definite  matrix  S has  the  Cholesky  decomposition 

S = L fT, 

where  L is  a lower  triangular  matrix.  Thus 


i 


Mz  = yLL^z 


L~1  Mz  = piTz  ' 

( 


PlTz 


m ••  rw 


(L"l  M L"'r)x  = px  , 


(23) 


where  x = l'^z  . (24) 

Therefore,  the  eigenvalues  of  M L“'^  are  pre£isely  the  eigenvalues  of 
equation  (19),  and  the  eigenvectors  x of  L”^  M and  the  eigenvectors 
z of  S“^M  are  related  by  equation  (24).  Note  also  that  equation  (23)  is  a 
Hermitian  eigenproblem,  since  L“^  M L~'^  is  a Hermitian  matrix.  It  is 
therefore  possible  to  solve  equation  (23)  by  numerical  methods  designed  for 
Hermitian  eigenprobl ems  which  explicitly  use  the  fact  that  the  eigenvalues 
are  real  (reference  2).  Therefore,  the  eigenvalues  computed  by  using  equa- 
tion (23)  will  always  be  real  as  required. 

This  computational  procedure  seems  to  require  a prohibitively  large  num- 
ber of  arithmetic  operations;  however,  reference  2 shows  that  the  computa- 
tions may  be  done  very  efficiently  because  of  the_special  structure  of  the 
matrices  involved.  For  example,  the  matrix  ML“^  can  be  computed 

(without  inverting  the  matrix  L)  by  using  only  2/3  n^  complex  multiplica- 
tions. This  compares  to  1/2  n^  complex  multiplications  in  the  computation 
of  alone  in  equation  (23).  In  terms  of  storage  required,  computation 
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time,  and  numerical  accuracy,  the  use  of  equations  (23)  and  (24)  is  prefer- 
able to  equation  (22). 

The  routine  in  reference  3 was  adapted  to  the  Hermitian  case,  using  the 
routines  available  in  reference  2 to  solve  the  eigenproblem  (23).  This 
routine  is  called  PENCLH,  and  its  listing  is  available  in  appendix  B.  The 
listings  of  the  routines  used  from  reference  2 are  not  available  because 
they  are  proprietary  information  under  terms  of  the  lease  arrangements  made 
with  International  Mathematical  and  Statistical  Libraries,  Incorporated. 
Finally,  it  is  pointed  out  that  the  routine  PENCLH  computes  all  the  eigen- 
values and  eigenvectors  of  equation  (/I),  and  not  merely  the  largest  eigen- 
value and  corresponding  eigenvector(s) . 


NUMERICAL  SOLUTION  OF  THE  EIGENPROBLEM:  INDIRECT  METHOD 


As  discussed  in  the  subsection  COMPUTER  STORAGE  PROBLEM,  the  drawback  to 
the  direct  method  is  excessive  computer  storage  for  large  arrays.  The  group 
coordinate  relaxation  (or  indirect)  method  overcomes  this  drawback,  but  at 
the  cost  of  computer  time  and  the  loss  of  ability  to  compute  the  lower  order 
eigenvalues/eigenvectors.  The  group  coordinate  relaxation  method  is  detailed 
in  reference  5 for  the  real  syirmetric  eigenproblem  Ax  = yx.  This  method  can 
be  extended  easily  to  the  Hermitian  eigenproblem 

Mz  = ySz  . (25) 


Although  the  method  can  be  extended  to  arbitrary  Hermitian  matrices  M and  S, 
wit^h^  S jsojsi t^ive^ definite k i.t  is  important  hei;e  to.  retain  the  structure  of  M 
and  S as  given  by  equations  (20)  and  (21).  The  reason  is  that  the  Hermitian 
forms  of  M and  S can  be  evaluated  directly  without  knowledge  of  any  of  the 
entries  in  the  matrices  M and  S.  This  is  the  fact  that  allows  the  computer 
storage  problem  to  be  overcome. 

The  following  notation  will  be  very  useful.  Define  the  basis  vectors 


ej  = <1  0 0...0  0>^  (26a) 

e^  = <0  1 0...0  0>^ 


• • 

f T 

e » <0  0 0. . .0  1>  . 
n 

^ • Note  that  each  of  these  vectors  is  of  length  n.  To  define  vectors  e^i  for 

m > n + 1,  we  first  set 


n , if  m is  an  integer  multiple  of  n (34b) 

m - [m/n]D,  if  not 

where ^ ^denotes  the  greatest  integer  function.  Since  equation  (26b)  requires 
1 < t < n,  we  can  now  define 


I 


t 


) 

t 


I 

■J 

1 

I 
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ra>n+l.  (26c) 

In  other  words,  we  have  defined 

“ ®2n+l  “ ••• 

®2  ° ®n+2  “ ®2n+2  “ 


I 


f 


f 


■2n 


= e 


3n 


Before  the  group  coordinate  relaxation  algorithm  can  begin,  two  items 
must  be  specified.  First,  an  initial  guess 


®(0) 


(0).T 
a > 
n 


(27) 


at  the  optimal  element  excitation  vector  is  required.  The  vector  a(o)  can- 
not consist  of  only  zero  entries,  but  it  is  otherwise  completely  arbitrary. 
Second,  it  is  decided  in  some  manner  to  work  with  subarrays  of  the  full  array 
of  size  r ^ 1.  It  will  be  shown  below  that  choosing  to  work  with  subarrays 
of  size  r will  mean  that  generalized  eigenproblems  of  size  r+1  will  have  to 
be  solved,  so  computer  storage  plays  an  important  role  in  the  choice  of  r. 
Another  important  consideration  is  computation  time.  In  general,  the  larger 
r is  taken  to  be,  the  faster  optimum  excitations  of  the  full  array  can  be 
computed . 

The  group  coordinate  relaxation  algorithm  is  most  easily  described  by 
exhibiting  the  first  two  steps  of  the  algrithm.  From  these  steps  it  is  easy 
to  see  the  general  procedure.  In  the  first  step,  we  seek  to 


-T 
X Mx 

maximize  

xeQp  x'^Sx 


(28) 


where  Qq  is  the  vector  space  of  dimension  r+1  whose  general  element,  x,  can 
be  written  in  the  form 


X 


'0  ^0) 


+ c 


1 1 


+ c e 


r r 


(29) 


for  some  complex  constants  c.,  Cp...,c  . It  is  shown  below  that  (28)  is  a 
ratio  of  Hermitian  forms  in  the  parameters  Icq,  Cp...,c^}.  Therefore,  by 
Theorem  1,  the  solution  of  (28)  requires  solving  an  eigenproblem  of  size  r+1. 
Let 


I 


*(1) 


*^0  ®(o)  ‘'i  ®i 


+ c e 
r r 


(30) 
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be  a vector  for  which  the  maximum  (28)  is  attained.  This  completes  the  first 
step.  In  the  second  step,  we  seek  to 


-T 

. . X Mx 

maximize  , 

-T 

xeQj  X Sx 

where  Qj  is  the  vector  space  of  dimension  r+1  whose  general  element,  x, 
be  written  in  the  form 

X = c_  a,,v  + c,  e , + ...  + c e„ 

0 (1)  1 r+1  r 2r 


(31) 


can 


(32) 


for  some  complex  constants  c^,  Cj^,...,c^.  Since  (31)  is  again  a ratio  of 
Hermitian  forms  in  the  parameters  {cq,  Cj,...,c^},  we  solve  an  eigenproblem 
of  size  r+1  to  compute  a vector 


^2)  " ‘'O  ^(1)  ^ '^l  ®r+l  ^ 


+ c e„ 
r 2r 


(33) 


for  which  the  maximum  (31)  is  attained.  This  completes  the  second  step. 
Continuing  in  this  fashion  defines  the  group  coordinate  relaxation  algorithm. 


We  see  that  this  algorithm  "cycles"  through  the  entire  array  using  sub- 
arrays  of  size  r.  This  is  because  the  basis  vectors  lej^}  are  defined  to 
"cycle"  regularly  through  the  vectors  {e^,  e2,...,en}.  Also,  if  r does 
not  divide  n evenly,  then  each  individual  element  belongs  tg  a.  ijumber^of  .dif- 
ferent subarrays  as  the  computation  proceeds.  In  other  words,  if  r does  not 
divide  n,  the  entire  array  is  not  subdivided  into  disjoint  subarrays. 

The  group  coordinate  relaxation  algorithm  generates  a sequence  of  vectors 
^(O)’  ^(1)»  ®(2)»’"’  which  converges  to  the  eigenvector  corresponding 
to  the  largest  eigenvalue  of  equation  (19).  Convergence  is  assured  regard- 
less of  the  starting  vector,  with  some  highly  unlikely  exceptions.  These 
exceptions  are  easy  to  state.  If  any  of  the  computed  vectors  {a(o), 

3(1),  a(2),--.}  is  precisely  an  eigenvector  of  equation  (19)  which  corre- 
sponds to  an  eigenvalue  that  is  not  the  largest  eigenvalue  of  the  equation, 
then  the  group  coordinate  relaxation  method  will  not  move  from  this  eigenvec- 
tor. Numerical  roundoff  error  will  probably  prevent  this  in  practice.  For 
further  discussion  and  for  a convergence  theorem  whose  proof  can  be  extended 
to  the  present  situation,  see  reference  5.  For  an  application  of  these  math- 
ematical methods  to  other  problems,  see  reference  8. 


An  important  feature  is  that  the  last  computed  vector  a(]^)  gives  a 
larger  DIBCFp  than  the  previous  vector  a(|^_i).  This  is  easy  to  see  by 
observing  the  ratios  (28)  and  (31). 


Another  very  useful  observation  is  that  the  algorithm  requires  knowledge 
of  only  a(^^)  to  compute  a(]^+i).  This  means  that  if  computation  must  be 
interrupted  for  any  reason,  it  is  necessary  to  store  only  the  last  computed 
vector  in  order  to  restart  computations. 
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We  conclude  this  section  by  an  examination  of  the  maximum  (28).  Every- 
thing that  is  said  of  (28)  is  easily  translated  to  the  maximum  (31),  as  well 
as  all  the  other  maxima  required  in  the  group  coordinate  relaxation  algo- 
rithm. Note,  first,  that  putting  (29)  into  (28)  gives  the  identity 


xGQ„  -T., 
^0  X Sx 


^ z^Bz 


T 

where  z=<c„,c,,  ...,  c>.  G=rg..]  and  B =rb..1are  (r+1)  x (r+1)  Hermi  t i ; 

0’  1 r •-  ir 

matrices  whose  general  entries  are  given  by 


®00  ^(0)”^(0) 


®0j=^0)”^j  ’ j = ^ " 


g,  . = e,  Me.  , k,  j 
kj  k J 


'’OO  ^(0)^  ^(0) 


% ^(0)  ^ ’ .1  " 


'’ko  ®k  ^ ^(0)  ’ 


% = % • 


Thus  the  entries  of  G and  B are  computable  from  the  Hermitian  forms  of  M and 
S,  respectively.  Let  Vq(9,(}>)  be  the  field  pattern  of  the  entire  array  for 
the  excitations  a(Q).  Then  we  have,  explicitly. 


(q,^)€'3  ni\ 


|N(0,(}))VQ(0,(t>)l 


I 


TR  5889 


’kO  = 


k = 1 , . . . , r 

(36) 

®0k  ~ ®fc0  ’ 

k = 1 , . . . , r 

(37) 

®kj  ’ \j  ’ 

k,  j = 1 , . . . , r 

( 38) 

(39) 


whore  is  given  by  (20),  and  similarly, 


kO  = 


'’ok  “ ^kO 


k=l , . . . , r 
k = I , . . . , r , 


= ®kj  ’ j = 1*  •••.  r . 

.wbftr<?tS(^i  i-s  -given  iiy  'f21).’  Be'cause  VQrO,^))  can  be  computed  easily  for 
each  (e,(^),  we  see  that  equations  (35)  through  (42)  can  be  computed  effi- 
ciently in  terms  of  time  and  core  storage  requirements.  Now,  by  using 
Theorem  1,  we  see  that  the  maximum  of 


(40) 

(41) 

(42) 


— T 
z Gz 

-T 
z Bz 


is  achieved  by  any  vector 


•?'  > 
r 


which  is  an  eigenvector  of  the  largest  eigenvalue  of  Gz  = nBz.  Thus,  from 
(34),  we  see  that 


®(l)  = ^0  3(0)  + ^1  «!  + •• 
is  a vector  for  which  the  maximum  (28)  is 


+ c. 


e 


r 


attained . 
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EXAMPLES 


EXAMPLE  1:  A COMPARISON  WITH  DOLPH-CHEBYSHEV  DESIGN 

This  example  serves  two  purposes.  First,  it  provides  a comparison  with 
the  Dolph-Chebyshev  line  array  design.  Second,  it  gives  some  insight  into 
the  nature  of  the  lower  order  eigenvalues/ eigenvec tors . 

Suppose  that  we  have  a line  array  of  15  elements  which  lies  along  the  y- 
axis  (see  figure  1)  with  equal  spacings  of  0.5  wavelength,  where  the  wave- 
length X=l.  (The  units  of  length  are  irrelevant.)  Thus,  if  the  first 
element  lies  at  the  origin  with  coordinates  (0.,  0.,  0.),  the  15-th  element 
has  the  coordinates  (0.,  7.,  0,).  It  is  well  known  that  any  line  array  has 
a field  pattern  with  cylindrical  symmetry  about  the  array  axis.  Therefore, 
the  most  effective  choices  for  and^  would  seem  to  be  cy  1 i ndr  i cal  1 y 
symmetric.  We  define"^  to  be  the  set  of  all  directions  that  lie  within  8 
degrees  of  a normal  to  the  y-axis,  and  we  define  A3  to  be  the  collection  of 
all  other  directions.  Hence  is  16  degrees  "wide."  The  ambient  noise 
field  is  assumed  to  be  flat,  and  the  individual  elements  are  assumed  to  be 
omnidirectional.  Finally,  considering  the  construction  of  the  Fuller  points 
■^p,  we  choose  p = 24. 

The  above  data  completely  define  the  DIBCFp  array  problem.  The  rest 
is  understanding  the  computer  program  usage  (see  appendix  A).  In  appendix 
B,  a listing  of  the  entire  program  required  for  this  example  is  given.  The 
resvil  ts  of  the  execution  are  given  in  table  5.  Also,  the  field  pattern  (in 
the  xy-plane)  is  given  in  figure  4. 

0 
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Figure  4.  Field  Patterns  for  Excitations  in  Table  5. 
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The  Dol ph-Chebyshev  excitations  are  designed  exclusively  for  half-wave-  »• 

length  equispaced  line  arrays  with  omnidirectional  elements.  For  a given 
number  of  elements,  the  Dol ph-Chebyshev  excitations  depend  only  upon  the 
steered  direction  and  on  the  specified  sidelobe  level.  For  a broadside 
(i.e.,  steered  normal  to  the  line  of  the  array)  15-element  array,  the  Dol ph- 
Chebyshev  excitations  for  a 28  dB  si  delobe- level  field  pattern  are  given  in 
table  5.  The  corresponding  field  pattern  is  shown  in  figure  4. 

We  note  that  the  mainlobe  shape  of  the  Dol ph-Chebyshev  array  and  the  \ 

DlBCFp  array  are  indistinguishable  in  figure  4.  The  only  difference  lies 

in  sidelobe  structure.  We  can  also  see  that  by  sacrificing  approximately  , 

3 dB  in  the  sidelobe  nearest  the  mainlobe,  all  the  remaining  sidelobes  can  I 

be  made  smaller  than  the  overall  28  dB  sidelobe  level  of  the  Dol ph-Chebys hev  J 

array . 

Table  5.  Excitations  for  15  Element  Equispaced  I 

Line  Array:  Dol ph-Chebyshev  Vs  DIBCF24  ^ 


Element  No. 

Dol ph-Chebyshev 

DIBCF24 

1 

.34371 

.25687 

2 

.35775 

.39520 

3 

. 50403 

.54024 

4 

.65338 

.68290 

5 

.79108 

.81242 

6 

.90242 

7 

.97487 

.97920 

8 

1.00000 

1.00000 

9 

.97487 

.97920 

10 

. 90242 

.91188 

11 

.79108 

.81242 

12 

.65338 

.68290 

13 

. 50403 

. 54024 

14 

.35775 

.39520 

15 

.34371 

.25687 

, What  about  the  lower  order  eigenvalues?  The  first  four  eigenvalues/ 

I eigenvectors  are  listed  in  table  6.  (Note  that  the  eigenvector  of  D j in 

! table  6 is  also  given  in  table  5,  but  is  normalized  differently.)  Also,  the 

corresponding  field  patterns  are  given  in  figures  5 through  8.  We  r-"sark 
only  that  the  field  pattern  for  the  largest  eigenvalue  P]  nas  no  nulls  in 
the  mainlobe  region'Wl^,  whereas  the  field  pattern  for  P2  null  in 

n-  two  nulls  for  pq,  and  3 nulls  for 
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Table  6.  The  Four  Largest  Eigenvalues/Eigenvectors  of  Example  1 


Element  No. 

Ml  = .9894 

M2  * .8206 

M3  - .3231 

M4  = .0383 

1 

.0916 

.2755 

.4486 

.5106 

2 

. 1409 

.3158 

.3666 

.2141 

3 

.1927 

.3289 

.2442 

-.0356 

4 

.2436 

.3112 

.1066 

-.2042 

5 

.2898 

.2675 

-.0242 

-.2664 

6 

.3252 

.1929 

-.1413 

-.2454 

7 

.3492 

.1030 

-.2097 

-.1391 

8 

.3567 

.0000 

-.2403 

.0000 

9 

.3492 

-.1030 

-.2097 

.1391 

10 

.5252 

-.1929 

-.1413 

.2454 

11 

.2898 

-.2675 

-.0242 

.2664 

12 

.2436 

-.3112 

.1066 

.2042 

13 

.1927 

-.3289 

.2442 

.0356 

14 

.1409 

-.3158 

.3666 

-.2141 

15 

.0916 

-.2755 

.4486 

-.5106 

EXAMPLE  2:  A 105-ELEMENT  CYLINDRICAL  ARRAY 


This  example  illustrates  the  use  of  subarrays  (i.e.,  the  group  coordi- 
nate relaxation  method)  for  computing  optimum  DIBC  with  limited  computer 
storage.  We  select  an  array  with  105  elements  arranged  around  a cylinder. 
Specifically,  we  first  construct  7 rings  of  15  elements  each,  and  then  place 
the  axis  of  each  of  these  rings  along  the  x-axis  (see  figure  9).  The  exact 
positions  (and  element  numbers)  are  given  in  table  7,  where  the  units  of 
length  are  such  that  the  wavelength  X = ]. 


z 

* 


RINGS  OF 
ELEMENTS 


Figure  9. 


Arrangemenl.  of  Elements  in  Example  2. 


Each  element  of  this  array  has  a hemispherical  field  pattern  defined  in 
the  following  manner.  We  conceive  of  the  array  as  being  supported  by  a 
(transparent)  cylinder.  Through  each  element,  we  pass  ^a'^plane  parallel  to 
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Table  7 

. Coordinates  of  Elements  in  Example 

2 

Element 

Coordinates 

No. 

> 

c 

y 

z 

1 

.0000 

.7642 

.0000 

2 

.3337 

.7642 

.0000 

3 

.6674 

.7642 

.0000 

4 

1.0011 

.7642 

.0000 

5 

1.3348 

. 7642 

.0000 

6 

1.6685 

.7642 

.0000 

/ 

2.0022 

. 7642 

.0000 

8-14 

As  above 

.6982 

.3108 

15-21 

.5114 

. 5679 

22-28 

.2362 

.7268 

29-35 

-.0799 

. 7600 

36-42 

-.3821 

.6618 

43-49 

-.6183 

.4492 

50-56 

-.7475 

.1589 

57-63 

-.7475 

-.1589 

64-70 

-.6183 

-.4492 

71-77 

-.3821 

-.6618 

78-84 

-.0799 

-.7600 

85-91 

.2362 

-.7268 

92-98 

.5114 

-.5679 

99-105 

.6982 

-.3108 

the  cylinder  axis.  The 

field  pattern  of  an 

element  has  unit 

response  on  the 

side  of  the  plane  which 

DOES  NOT  contain  the 

cylinder  and  has 

zero  response 

on  the  side  that  DOES  contain  the  cylinder. 

We  assume  that  the  ambient  noise  field  is  flat.  Also,  we  choose  p = 32 
in  the  definition  of  the  Fuller  points  «^p. 

The  mainlobe  region*^  is  defined  as  a "half  cone"  lying  above  the  posi- 
tive x-axis.  Specifically,  consider  the  solid  cone  with  axis  lying  along 
the  positive  x-axis,  with  its  vertex  at  the  origin,  and  with  a vertex  angle 
of  40  degrees.  The  xy-plane  slices  this  cone  into  two  equal  parts,  and  the 
mainlobe  region*^  is  defined  to  be  that  part  of  the  cone  which  lies  above 
(i.e.,  points  having  positive  z-coordinates ) the  xy-plane.  The  sidelobe 
region /tf  is  defined  to  be  the  set  of  all  directions  that  are  not  in  the  main- 
lobe region  There  is  no  ignored  regional  in  this  example. 

With  all  the  above  choices,  the  DIBCFp  array  problem  is  completely 
specified.  To  solve  the  entire  problem  at  once  on  the  UNIVAC  1108  requires 
more  core  storage  than  is  available  without  resorting  to  indirect  address- 
ing or  some  other  scheme,  so  we  use  subarrays  to  maximize  the  DIBCFp  of 
the  entire  array.  It  seems  best  to  use  as  many  elements  as  can  be  easily 
handled  with  the  computer  storage  available,  so  in  this  case  we  choose  69 
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I 

\ elements;  that  is,  roughly  2/3  of  the  array.  This  turns  out  to  require  only 

about  45,000  words  of  computer  storage.  To  effect  this  computation,  using 
the  program  in  appendix  B,  requires  that:  (1)  we  use  the  computer  program 
as  if  we  were  going  to  solve  the  entire  problem  without  group  coordinate 
relaxation,  and  (2)  we  add  only  2 additional  cards,  namely,  NRELAX  = 69+1 
and  EPSI  = .001. 

The  group  coordinate  relaxation  scheme  required  roughly  1650  seconds  per 
iteration,  and  5 iterations  in  all.  Thus  total  computation  time  was  roughly 
2 1/4  hours.  Table  8 gives  the  final  (optimal)  set  of  element  excitations. 
Figure  10  gives  the  vertical  field  pattern,  and  figure  11  gives  the  horizon- 
tal field  pattern,  for  these  excitations.  We  point  out  that  in  the  field 
patterns  in  these  two  figures  have  abrupt  jumps  because  the  individual  ele- 
ment field  patterns  also  have  sharp  jumps  due  to  their  assumed  hemispherical 
field  patterns.  (These  field  patterns  were  computed  by  the  field  pattern 
program  described  in  reference  9.) 


Table  8.  Optimum  Excitations  for  Example  2 


I 

N«). 

M.if»n  i f iiH*' 

Ph.isp 

K 1 *'mrnt 

No. 

Mij»n  if  i\<1c 

Phasp 

I 

.01 497 

-7.04140 

44 

.049! 7 

- .20842 

2 

. 04  4 44 

1.  444  78 

44 

.04)07 

-3.12147 

.07401 

-1.77318 

44 

.01413 

.24034 

u 

.091  7^1 

7.1474/* 

47 

.01474 

-7.08431 

''i 

.08 ->04 

- .72447 

48 

.04234 

1 . 42329 

6, 

.04214 

2.48449 

49 

.09072 

- 1 .29944 

7 

.07174 

- .14344 

40 

.10997 

2.13497 

R 

- . 7?84/> 

4; 

.09774 

- . 74037 

•) 

.OISHH 

2.44704 

47 

.041  '*2 

2.47418 

10 

.04133 

- .304«8 

43 

.07414 

- .14489 

1 ] 

.0734 

•>.914  74 

44 

.01874 

-2.47888 

P 

.04403 

.000  79 

44 

.04833 

.97472 

n 

.0380  ) 

-2.97334 

44 

.10412 

-1.80774 

.01308 

..*4931 

47 

. 13777 

1.44409 

.00348 

- .013'’3 

48 

.12180 

-1.  18747 

1^ 

.017  78 

3. 1 3121 

49 

.07979 

’.7483? 

l 7 

.07704 

.29844 

70 

.0321? 

- .48940 

iH 

.07717 

-’.40443 

71 

.03444 

2.  70718 

10 

.07310 

. 987/*'4 

77 

.08437 

- .0441? 

.01400 

- 1.  /I-JV. 

73 

.12944 

-2.8481? 

.00470 

1 . 78944 

74 

.14744 

.41 788 

?? 

. 01037 

1.48  324 

74 

. 1 1 737 

-’.1  7831 

21 

.07424 

-1.20497 

74 

.noo-*? 

1.34010 

*>4 

.04414 

2.09737 

77 

.024  74 

- 1.30499 

2^ 

.04783 

- .90479 

78 

.07023 

2.82734 

24 

.04388 

7.37497 

79 

.1884’ 

.01  770 

“>7 

.03470 

- .01774 

80 

. 3089? 

-2.83428 

2H 

.01410 

2.  70344 

81 

.34744 

.49043 

20 

.03.’'»7 

I . 3'>444 

82 

. 30448 

-2.24194 

30 

.08731 

- 1 . 4r)<)9ft 

83 

.18437 

1.18412 

31 

.14424! 

1 . 741  47 

84 

.04783 

-1 . 494  34 

3? 

. 1 7443 

- 1 .’0741 

84 

.04400 

2.42014 

33 

. 1 4 74  3 

2. 17444 

84 

. 1 ’949 

- .31104 

3'-i 

. 101  04 

- .'illW 

87 

.71923 

3. 11424 

3'i 

.04090 

2.47482 

88 

.24487 

.244  34 

34 

.00107 

1 . 8 0)0  1 

89 

.70337 

-2.  49481 

3 7 

.00307 

-!.4790> 

9f» 

. 11910 

.84311 

3*3 

. 00847 

1 . 48840 

01 

.04179 

-I.H444R 

30 

.0! 444 

- 1 . 30O->  \ 

()•* 

.02424 

2.88974 

40 

.01414 

1.04430 

93 

.04438 

.72940 

41 

.01282 

- .04279 

'•4 

.10232 

-2.41417 

4? 

.00443 

2.49084 

94 

.11889 

1.00483 

4 3 

.00811 

. 10499 

04 

.10393 

-1,74944 

44 

.0730) 

-2.94378 

97 

.04423 

1.74744 

44 

.04148 

.42414 

98 

.02474 

- .94427 

40 

.04011 

‘7.47388 

99 

.01431 

-:.  44878 

47 

.0/.477 

.02478 

100 

.044f*4 

1.14241 

4« 

.077»4 

- t.O'’483 

loi 

.081 78 

-1.43309 

40 

.01048 

1.40198 

10? 

. 10300 

1.87940 

. 01 744 

-1.14074 

103 

.09444 

-1,01412 

41 

.03073 

7.24479 

104 

.041’4 

2.42174 

42 

.047V« 

- .44477 

104 

.0'»421 

- .41890 

43 

.07080 

2.  70701 

I- 
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Figure  11.  Horizontal  Beam 
Pattern  for  Example  2 with 
Excitations  Given  in  Table  8. 


ANGLE  MEASURED  FROM  y-AXIS  (deg) 


I 

I 
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This  method  creates  a steadily  increasing  sequence  of  estimates  for  the 
largest  eigenvalue.  Since  there  were  5 iterations,  there  were  5 estimates,  anvi 
these  are  given  in  table  9.  Based  on  the  plots  in  figures  10  and  11,  it  would 
seem  that  the  field  pattern  cannot  be  improved  significantly.  In  other  words, 
even  if  we  had  chosen  the  variable  EPSI  to  be  even  smaller,  we  would  have 
forced  additional  (expensive)  iterations  of  the  algorithm  which  would  probably 
not  have  improved  the  field  patterns. 


Table  9.  Group  Coordinate  Relaxation  Estimates 
of  Largest  Eigenvalue  for  Example  2 


T terati on  No. 

Estimate  of  Largest  Eigenvalue 

1 

.91427 

2 

.96100 

3 

.96374 

4 

.96532 

5 

.96559 

EXAMPLE  3:  EFFECTS  OF  SAMPLING 

The  first  two  examples  did  not  mention  the  effects  of  sampling  on  the  field 
patterns.  Specifically,  the  parameter  p in  the  definition  of  the  Fuller  points 
U,p  determines  how  finely  we  have  sampled  all  spatial  directions.  Hence,  the 
parameter  p has  a strong  influence  on  the  resulting  field  patterns.  In  partic- 
ular, if  p is  not  sufficiently  large,  then  it  is  possible  for  the  optimal 
(DIBCFp)  field  pattern  to  have  a split  beam. 

We  illustrate  this  effect  by  systematically  varying  p in  the  array  of  exam- 
ple 1,  but  for  a different  choice  of7^^  and /</ . Here,  we  define  7*1  to  be  the 
collection  of  all  directions  whose  projection  onto  the  xz-plane  lies  within  ^8 
degrees  of  the  z-axis.  Specifically,  the  direction  corresponding  to  direction 
cosines  (a,P,Y)  lies  infri  only  if  |a/<f afi  + ' 1 £ sin  8°.  In  other 

words , “Hj cons  is ts  of  all  directions  contained  "between"  the  two  planes  inter- 
secting the  yz-plane  at  the  angles  of  +8  degrees  and  -8  degrees.  The  sidelobe 
region;0  consists  of  all  remaining  directions,  so  there  is  no  ignored  regions!. 
Optimal  excitations  for  several  choices  of  p are  given  in  table  10.  The  field 
patterns  for  p = 24  and  p = 16  are  given  in  figure  12.  We  do  not  present  the 
field  patterns  for  p = 32  and  p = 40,  because  they  are  so  similar  to  p = 24. 

Table  10.  Effects  of  Sampling  on  Excitations  for  Example  3 


R1 

.... 



^2 

- * ' 

No. 

1ft 

24 

32 

40 

1 

.2933 

.1273 

.1322 

.1177 

2 

-.00ft4 

.1676 

. 1 792 

.1662 

3 

.384ft 

.2207 

.2406 

.2232 

4 

-.1043 

.2497 

.2613 

.2381 

.2990 

.2884 

.2946 

.2934 

ft 

-.2830 

.3067 

.2991 

.3099 

7 

.1733 

. 3337 

.3134 

.3239 

8 

-.3280 

.3346 

.3117 

. 3279 

9 

.1733 

.3337 

.3134 

.3239 

10 

-.2830 

.3067 

.2991 

.3099 

II 

.2990 

.2  884 

.2946 

.2934 

12 

-.1043 

.2497 

.2613 

.2381 

n 

.3846 

.2207 

.2406 

. 2232 

14 

-.0064 

.1676 

.1792 

.1662 

13 

.2933 

. 1273 

.1322 

. 11  77 

Lar||»<it 

ci  grnval  ua  • 

M .1149 

.1243 

.1163 

.1223 

No  . of  Fii  11  vr 

pdi  nt  A ,9  p 

2362 

3762 

10242 

16002 
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Figure  12.  Field  Pattern  for  Example  3 for  p = 16. 


Table  10  also  shows  the  effects  of  oversampling.  Note  that  the  optimal  [ 

excitations  for  p > 24  are  all  similar,  but  they  do  not  seem  to  be  converg-  1 

ing  to  an  optimal  set.  This  is  probably  due  to  the  buildup  of  numerical  } 

round-off  error  in  the  required  sums  (i.e.,  equations  (20)  and  (21)),  but  it  ! 

could  also  be  that  p must  be  chosen  even  larger  than  40  before  the  optimal 
excitations  give  the  appearance  of  convergence.  In  any  event,  the  import- 
ance of  sampling  sufficiently  finely  is  clear,  but  evidently  oversampling 

wastes  time  and  increases  the  numerical  round-off  error  in  the  computed  opti-  t 

mal  excitations. 


EXAMPLE  4:  TIME  AND  ACCURACY  IN  THE  INDIRECT  METHOD  , 

f 

It  is  clear  from  the  definition  of  the  indirect,  or  group  coordinate  ' 

relaxation  method,  that  the  size  of  the  subarrays  used  and  the  stopping  cri-  ^ f 

teria  for  the  iteration  procedure  both  have  significant  effects  on  numerical  j 

accuracy  of  the  computed  excitations  and  on  the  time  required  to  compute 
them.  This  example^ il lustrates  how  numerical  accuracy  and  computation  time 
depend  on  both  these  parameters. 

We  consider  a line  array  of  25  elements  that  lies  along  the  y-axis  with  . 

equal  spacings  of  0.5  wavelength,  where  the  wavelength  X • 1.  Thus  the  co-  ■ i 

ordinates  of  the  first  and  last  elements  are  (0.,  0.,  0.)  and  (O.,  12.,  0.)  I 

respectively.  We  select  the  mainlobe  region  >L  to  be  the  set  of  all  direc-  • 

tions  that  lie  within  5 degrees  of  a normal  to  the  y-axis,  and  we  define 
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to  be  the  set  of  all  other  directions.  There  is  no  ignored  region  . The 
ambient  noise  field  is  flat  and  the  individual  elements  are  assumed  omnidi- 
rectional. Finally,  we  select  the  Fuller  points3*2g*  All  the  above  com- 
pletely define  our  problem. 

Table  11  shows  the  number  of  iterations  required  for  various  choices  of 
NRE1.AX  (i.e.,  subarray  size)  and  EPSI  (i.e.,  stopping  criteria).  See  appen- 
dix A.  As  can  be  expected,  the  number  of  iterations  required  increases wi th 
decreasing  EPSI  and  decreases  with  increasing  NRELAX.  Also,  the  computation 
time  per  iteration  increases  with  NRELAX. 


Table  11.  Number  of  Iterations  Required  in  Example  4 


NRELAX- 1 

10-3 

EPSI 

10-A 

10-5 

Tirae/I terat ion  (sec) 

5 

9 

13 

24 

83 

10 

5 

10 

13 

95 

15 

3 

4 

6 

115 

20 

3 

3 

3 

140 

25 

1 

1 

1 

172 

An  important  concern  is  the  numerical  accuracy  of  the  computed  excita- 
tions. This  is  particularly  important  in  light  of  the  fact  that  numerical 
computation  of  eigenvectors  by  any  method  is  considerably  less  stable  than 
the  numerical  computation  of  eigenvalues.  Table  12  shows  the  results  ob- 
tained for  NRELAX- 1 =’  5 by  stopping  after  the  first  four  complete  passes 
through  the  array,  i.e.,  for  iterations  5,  10,  15,  and  20,  respectively. 

The  exact  results  are  included  also.  The  field  patterns  corresponding  to 
excitations  of  iteration  5 and  the  exact  excitations  are  shown  in  figure 
15.  Note  that  at  the  end  of  iteration  5 the  field  pattern  already  possesses 
si  delobes  in  the  correct  positions,  although  they  are  about  3 dB  higher  than 
in  the  field  pattern  of  the  exact  excitations.  Thus,  the  effect  of  later 
iterations  is  to  beat  down  the  sidelobes  while  maintaining  the  mainlobe 
beamwidth. 
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Table  12. 

Example  4 with 

Subarrays 

of 

5 Elements 

Element  No. 

Iteration  5 

Iteration  10 

I teration 

15 

Iteration  20 

Exact 

1 

.632 

.410 

.623 

. 793 

1.000 

2 

.813 

.622 

.887 

1.090 

1.339 

3 

.993 

.860 

1.173 

1.406 

1.692 

4 

1.183 

1.117 

1.428 

1.736 

2.057 

5 

1.394 

1.398 

1.800 

2.082 

2.436 

6 

1.446 

1.788 

2.218 

2.461 

2.  792 

7 

1.648 

2.079 

2.539 

2.795 

3.145 

8 

1.840 

2.355 

2.832 

3.095 

3.456 

9 

2.023 

2.611 

3.100 

3.366 

3.  732 

10 

2.  188 

. 2.842 

3.  328 

d.  592 

3.955 

11 

2.  541 

3.143 

3.  522 

3.783 

4.120 

12 

2.664 

3.298 

3.654 

3.904 

4.223 

13 

2.  749 

3.392 

3.722 

3.956 

4.254 

14 

2.806 

3.439 

3.737 

3.951 

4.223 

15 

2.809 

3.421 

3.677 

3.877 

4.120 

16 

3.011 

3.271 

3.559 

3.729 

3.955 

17 

2.951 

3.145 

3.396 

3.  541 

3.  732 

18 

2.834 

2.967 

3.  174 

3.298 

3.456 

19 

2.696 

2.  756 

2.929 

3.022 

3.145 

20 

2.489 

2.494 

2.631 

2.  700 

2.  792 

21 

2.012 

2.190 

2.290 

2.353 

2.436 

22 

1.762 

1.887 

1.957 

2.000 

2.057 

23 

1.511 

1.567 

1.632 

1.658 

1.692 

24 

1.260 

1.294 

1.313 

1.324 

1.339 

25 

1.000 

1.000 

1.000 

1.000 

1.000 

D max 

.9798313 

.9857357 

.9880019 

.9886998 

.9889890 

SUMMARY 


The  concept  of  Directivity  Index  with  Beamwidth  Control  (DIBC)  has  been 
defined  as  the  ratio  of  power  in  the  mainlobe  region  to  the  total  power  in 
both  the  mainlobe  and  the  sidelobe  regions.  A mathematically  and  numerical- 
ly tractable  method  for  the  computation  of  optimim  element  excitations  (i.e., 
excitations  which  maximize  DIBC)  is  presented.  A technique  known  as  group 
coordinate  relaxation  is  shown  to  be  an  effective  means  of  computing  optimum 
element  excitations  for  arrays  of  arbitrary  numbers  of  elements,  yet  it 
requires  only  nominal  core  storage.  Conceptually,  the  group  coordinate 
relaxation  technique  employs  subarrays  of  the  full  array  in  a systematic 
manner  to  optimize  excitations  of  the  full  array.  Four  examples  have  been 
included,  one  of  which  demonstrates  the  effectiveness  of  group  coordinate 
relaxation  for  a cylindrical  array  of  105  elements. 
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APPENDIX  A 

PROGRAM  USAGE 


The  user  is  required  to  write  between  2 and  4 routines  which  are  to  be 
used  Ln  conjunction  with  the  elements  in  the  program  file  STREIT*RELAX.  The 
setup  and  use  of  each  of  these  routines  is  easily  seen  from  the  listings  in 
appendix  b. 

SUBROUTINE  PHONE  - Given  a set  of  direction  cosines  (Ct,  B,  Y),  compute  the 
directional  response  of  every  individual  element  in  the  direction  (Ot,  8,  Y). 
If  all  elements  are  omnidirectional,  this  subroutine  need  not  be  written  by 
t he  us  er . 

FUNCTION  DNOISE  - Given  a set  of  direction  cosines  (Ot,  B,  Y),  compute  the 
magnitude  of  the  ambient  noise  field  in  the  direction  (Qt,  B,  Y).  If  the 
ambient  noise  field  is  flat,  this  function  subroutine  need  not  be  written  by 
t he  us  er . 


SUBROUTINE  REGION  - Given  a set  of  direction  cosines  (Ot,  B,  Y),  decide  (by 
user  prerogative)  if  this  direction  lies  in  the  mainlobe  region'Aj,  the  side- 
lobe  region  or  the  ignored  region^.  This  subroutine  must  always  be 
written  by  the  user. 


MAIN  PROGRAM  - This  routine  sets  up  the  array  problem  and  handles  all  the 
required  computer  allocations.  Because  of  the  computer  storage  problem,  all 
the  dimensional  storage  arrays  (not  to  be  confused  with  the  antenna  array) 
are  defined  by  means  of  one  PARAMETER  statement.  Therefore,  the  main 
routine  must  be  recompiled  for  each  different  array.  The  following 
variables  must  be  specified  within  the  main  routine: 


N 


Total  number  of  elements  in  the  antenna  array.  N MUST  BE 
DEFINED  BY  MEANS  OF  A PARAMETER  STATEMENT. 


NRELAX 


X(I) 

Y(T) 

Z(  I) 

WAVLEN 


I NOISE 


If  group  coordinate  relaxation  is  not  to  be  used,  NRELAX  = 
N.  If  group  coordinate  relaxation  is  to  be  used,  NREIAX  is 
1 plus  the  number  of.  elements  to  be  relaxed.  In  either 
case,  NRELAX  MUST  BE  DEFINED  BY  MEANS  OF  A PARAMETER  STATE- 
MENT. 


The  coordinates  of  the  I-th  element,  I = I,  2,  ...,  N. 


Wavelength  of  the  design  frequency  in  units  compatible  with 
the  element  coordinates. 

0 if  omnidirectional  noise  field 

1 if  noise  field  is  specified  by  subroutine  DNOISE. 
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IPHONE 


INDEX 


ALP HAD 

BETAD 

GAMMAD 

I STEER 


ALPHAS 
B ETAS 
GAMMAS 

NDIV 

RATIO 


COEF 


ICOEF 


0 if  every  element  in  the  antenna  array  has 
tional  field  pattern 

1 if  at  least  one  element  has  a directional 
specified  by  subroutine  PHONE. 


an  omnidirec- 
field  pattern 


I if  ordinary  directivity  index  is  to  be  maximized. 
0 if  not. 


The  (oi,  B,  y)  direction  cosines  of  the  direction  in  which  to 
maximize  DI,  provided  INDEX  =1.  If  INDEX  = 0,  these 
variables  are  not  needed. 


0 if  the  array  is  unsteered 

1 if  the  array  is  to  be  steered 


The  (OL , B,  Y)  direction  cosines  of  the  direction  in  which  to 
steer  the  array,  if  ISTEER  = 1.  If  ISTEER  = 0,  the  default 
values  are  ALPHAS  = BETAS  = GAMMAS  = 0. 


The  parameter  of  the  Fuller  points  .^p;  that  is,  NDIV  = p. 

On  return,  an  estimate  of  the  mainlobe  power  to  sidelobe 
power.  If  INDEX  = 1,  an  estimate  of  DI.  (This  ratio  can  be 
unre liable . ) 


On  input,  COEF  is  unimportant  if  NRELAX  = N.  If  NRELAX  < N, 
then  COEF  contains  the  initial  guess  at  all  N of  the  final 
excitations.  On  return,  COEF  contains  the  optimal 
coe  f f ic i ents . 


= is  the  number  of  coefficient  sets  to  print;  i.e.,  print  the 
ICOEF  largest  eigenvalue/eigenvector  sets.  1 ^ ICOEF  ^ N. 
(Cannot  be  used  for  NRELAX  ^ N.) 


If  NREFAX  < N,  then  the  following  additional  variables  are  needed: 


LOOPMX  - maximum  number  of  relaxation  iterations  to  perform. 


EPSI  * stopping  criteria.  Stop  if 


IPRINT 


I PUNCH 


old  eigenvalue  estimate 
new  eigenvalue  estimate 

0,  if  intermediate  printout  not  desired 
Ij  if  intermediate  printout  is  desired 

0,  if  no  card  output  is  desired. 

1,  if  all  newly  computed  excitations,  real  and  imaginary 
parts,  are  to  be  punched  in  FORMAT  ( 3( 14 , 2E1 1 . 6 ) ) . 


< EPSI. 


The  listing  of  the  program  in  appendix  B is  also  the  listing  of  the  code 
required  for  the  particular  DIBCFp  array  problem  of  example  1.  The  run 
stream  is 


(3  RUN  (user's  RUN  card) 

0 HDG  FIRST  EXAMPLE 
0 ASG,A  STRE IT* RELAX. 

0 ASG,A  IMSL*MATHLIB. 

0 FOR, IS  STREIT*RELAX.MAIN 
0 FOR, IS  STRE IT*RELAX. REGION 
0 FOR, IS  STREIT*RELAX. PHONE 
0 FOR, IS  STREIT*RELAX,DNOISE 
0 MAP, IS 

IN  STREIT  *RELAX. 

LIB  IMSL  *MATHLIB. 

END 
0 XQT 


NOT  NEEDED  IN  EXAMPLE  1 


A-3/A-4 
Reverse  Blank 
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APPENDIX  B 
PROGRAM  LISTING 


brHLlT»HE.LAX(l) 

1 t 

2 (. 

3 C 

<*  C 

6 C 

t>  C 

^ c 

ts  t 

’ c 

10  t 

11  c 

12  C 

U (. 

1«  (. 

lb  c 

lb  t 

17  (. 

1«  C 

1*^  c 

20  c 

21  l 

22  C 

23  C 

2“  L 

2b  C 

26  C 

27  (. 

26  C 

29 

30  (, 

31 

32 

33 
3m 
3b 

36 

37 

3«  c. 

39  c 

MO  ^ . 

41  c 

4^  C 

43 

^4 

4b 

4h 

47 

4H 

49 

bO 

bl 

bi 

b3 

b4 

bb 

bb 

b7 

bb 

b9 

bP 

bl 

b2 

b3  C 

b4  t 

bb  C 

bb  C 

b7  L 

bb 

o9 

7P 


• MAIN 

b'lBROUTlNE  main  wRlTTfN  ROY  L#  ST^eij,  NAVAL  UNDERWATrR  SYSTEMS 
center*  new  LONDON  LABORATORY#  NEw  LONDON,  rT  10  AUGUST  1977 

IHIS  main  program  sets  l-P  AN  ARRAY  CONP IGURa TTON  AND  l»SFS  GROUP 
COORDINATE  relaxation  TQ  SOLVE  POR  OPTI'Hjm  rLPMENT  CURRENTS. 

N = total  number  of  elements  in  the  array 

NRElAX  = 1 ♦ the  N'UMHER  of  element  CURHEmt^  to  •relax*  Im  faCM  step 

NOIV  S parameter  controlling  the  number  of  points  chosen  cm  the  sphere 
X(I)  2 X coordinate  of  I-TH  element 

Y(I)  s Y coordinate  OE  I-TH  ELEmfnT 

7(1)  = I coordinate  of  l-TM  ELEMENT 

wAVLEr  X i*<AVEL£UGTH  OF  DESIGN  FRFOUFNCY 

LOOf’Mx  s MAX  NUmbfR  of  RELAXATIONS  TO  PERFORM,  qn  RETURN,  OONr , 

EPSI  s STOPPING  CRITERIO* ...STOP  IF  TMWHOVrMPNT  IS  NOT  AT  LEAST  FPSI 
COEF  r INITIAL  GUESS  AT  flEVENT  CURRENTS,  RETURN,  FINAL  AMStiER, 
ratio  = ON  RFTUPN,  AN  estimate  Op  maINLORF  POWER  TO  SIPELPBE  POWFR 
If-Oisr  = 1 IF  NOISE  FIEL''  present#  = *’  IF  NOT 
TPMOUr  X I IF  DlprCTlONAL  ELF'^ENTS  ARE  USF^,  x 0 IF  uqT 
tNDLX  X 1 IF  directivity  InDFX  IC  TO  Bt  mavI^TZEP,  x 0 tF  NOT 
(ALPUAJ#HETAP#Gam-'AO)x  direction  in  which  ^0  maxi  *IZE  ■“•T  TF  INOFX  X 1 
rSTtEP  = 1 IF  Array  IS  STpEREO#  s n JF  mot 

(ALf>HAS#METAS#GAM*AS)  i -IRECTION  COSI  TS  /^F  ARRAY  IF  ISTFFR  x l 
IPRINT  X 1 rr  Ir.T«’R-1EDlATr  p»INTOUT  TS  pFSfRrQ,  r 0 JF  NOT 
TPU\Cm  X J IF  iNTERMEOlATf  COEFFICIFNTS  AOr  TO  PE  PUNCHEo,  r 0 TF  NOT 
ICOtF  X OF  COEF  SFTS  TO  PPTNT  C'NLY  fop  CASE  NzNRFLAX  > 

TmE  RfMAlNOEP  OF  THF  arrays  ARE  ifORK  A'^FAS 

PAf-A  iLTtfJ  uxlS»NRELAXxlb 

NAf AnETLK  N3x3*NRELAx 

tOf  Moi,  IkoISE#  IPuO^'E,  Avl^N,  I'iOEX,  AL^MAl»,nFTAO#GAVA0»  ICOfF. 

1 IbTEER»  ALPHYS,^JFTAC,f,A“MAS,NOIV#LOOP  •X#FPSI , IPRINT,  T^'U^'CH 

CO^  PlE  » AMAlN(NKELAX,’JRtLAX) , ASiDF ( NPFLA X » NoFL AX ) »COEF(N) , 
i VX (NR£LAX»NrFLaX) .:iRECTcN) ,oc (NpFi  AV) 

uT  *E  #Sl  ON  X(N)  # Y(N)  #7(N)  ,''L(fJ)  » WK  (*n)  #OCNPFi  AVI  , iNOX  (“PFLAX)  , 

I V(NRELAX> 


nluT.  CU.^E  FOO  SPECIFIC  A^pay  PROBLEM 


UO  3 Krl #N 

A(K>xr, 

Y(n  ) = )/2, 

\ tu  )xO, 

hAVLLNrl, 
iNf  Exx** 
ISTEERxO 
I'lribEzO 
iPnO.Jt  xO 
APR! #T=0 
iPUNCUrO 
tPbIxl .fc-J 
ALPHASxO, 
bETASS''. 

OA  vv.aSxO  • 

ALPHADxALPHAS 

oFTAUSGtTAS 

GAmmaOxoAMmAS 

KCEFxib 

fiOTVx?ii 


tNL  COrt  FOR  SPECIFIC  ARRAY  PROBlf** 


call  HCuAX(N,NRELAV#X#Y,7,C0EF#RATl0,aM^p•,,«JTDE,nL#0•V»VX»•X» 

1 iNnx#rirECT»['C) 

LUO 


•Hr<T#b  '•«££/♦* 


B-l 


XBlfi  PlfflE  IS  BE5SX  QUALITY  PRACIlCABa 

fBoa  ocrx  wraished  to  ddo  — ^ 
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1 

i 

i 


I 

f 


bTRt  IT«r4ELAX(l  ) .MLLAX 

1 C SURROUTINP  RELAX  A.RCTTFN  RY  ROY  L,  ST^EIT,  ^'AVAU  UNnfRWATrQ  SYSTEVS 

i C LFMER*  (JEv^  LONDON  LAHORAToRY,  NEw  LONDON,  rT  10  AUONST  1077 

3 C 

4 C SliRROUTlNE  RELAX  USES  RROl'P  COORDINATE  RELAyATIOM  TO  COMPUTE 

b C oPTiMUr*  ELEMENT  CURRENTS  . 

R c 

7 C N = NiiMbER  OR  ELEMENTS  IN  THE  ApoaY 

R C NWELAx  z \ ♦ THE  MUMOEH  OF  ELEMENT  C'lRRENTt;  to  »RFLAX» 

0 c Yd)  = A COORDINATE  OF  I-TH  ELEMENT 

K C Y(l)  r Y COORDINATE  OF  I-TH  ELEMENT 

11  C 7(1)  = ^ COOROIrjATE  OF  I-TH  ELEmemT 

12  C CCEE  z initial  GUESS  AT  tlEMEMT  rURRRNTS 

10  c PATIO  = ON  RETUPhi,  AN  fSTI  ’ATE  OF  MfiTNUORE  '^owEP  TO  SIDPLOME  pc*(ER 

14  L AVAIU  = -MATRIX  OF  QUADRATIC  RORm  rOR  MAJNL^RP  REGION 

l‘j  c ASIDE  = MATRIX  OF  QUADRATIC  FORM  COP  SlOEL-^'^F  RFGIOrj 

Ih  L TmE  Hr-^PlNOEP  OF  THC  ARRAYS  ARE  V'ORX  A^^fAS 

17  C 


Irt 

IN 

2'' 

21 

tZ 

23 

24 

26 

27 

26 

30 

31 

32 

33 

34 
3b 
3b 
37 
3b 
39 

41 

4<; 

h3 


c 


G"*^R0UT  I if  RELAX(N,NRrLAX,  v,Y,Z»rOEP»RAriO,  .‘^AIN,ASlDE,rL»'',V»VX, 
I wx#  rJDX»''I9ECT*nr) 


->! VE  JSIO.«  X(  1 ) ,Y{  n ,7(1 ) ,:  J.(l)  ,p(^,  12)  » '(3#i2)  »D(l)  tr'OXd  ) ,WK  (1 ) , 
i V(  I ) »NTI  ‘E  (2) 

tovpLty  A ■«aUj(NREL  aX,  I ) , A«IPE  (^IRF|  ax,  1 ) »VX(*  relax,  1 ) ,rOEE(  1 ) , 

1 tIRECT(l),oC(l) 

LOvvor,  INOlSEi  IPhO*  E ,wA\Lf  N,  If  DFX,  ALPHA  ’,PC-»AO,GAMHAl),  ICO^F  , 

1 rjTEE«,ALPHAS,PETAf  ,GA‘-*MAS,NOTV,L00P*-X»ERsI  ,IPWI»IT,  I»n'^  Cm 

LO*  •^0'i/LATRA/ni2LA--,  ITIfiDT 
uATA  ((K(I,J)»l=l,^),jSl,12)/0,,n^,l,, 


1 
C 

3 

H 

G 

b 

f 

u 

'•* 

1 

2 


,804427191  ,0. » ,4*i7213S9S» 

,276393^02#  ,'*506nOR0fl,  ,44  721  3bP*>, 

,7?36«f.7RA»  •62^7^1 1 l2,-,44  7215Sqs, 

, 72360^798, -,62b^31 1 12,-,44721 

•P*,0»— 1,# 

-,7236rf79«»  .b2S«*31  112»,447213^oS, 
-,723f.bG7OR,-,c2«77iii2,,4472l3E0S, 
,2783Q'2'*2,-.flp'>x.S''R0A,  ,4472ns9b. 

-, 276307202 » ,86O«.Sf'80P,-, 44721 'SON, 
-,A94427101 , ,0,-,447213S9S» 

-,  27630  320  2 *-,P*'*S70  80P,-,  4 4 721  3b”6/ 


RPIriT 

^ORMAT(lMl,•  Tut  element  P'^SITIONS  IN  C'RTErIA^r  COORDINATES  A^E  *,• 
1/,  lr>X.  »tLT,  fin.*,  7X,  fX  • , 16X,  *y*  , ISX* ’Z*  ) 

RWINT  219.  (K,  X (k  ) , Y(t< ) ,7  (X  t ,*<rl  *»•) 

219  E0(  '1AT(  i2P.3El6,M) 


4b 

4b 

‘♦T 

40 

49 

bO 

bl 

b2 

b3 

b4 

bb 

b6 


PRINT  i2D.»AvLEN 

?..b  fOfvUAK/,*  wfAVELrnr.TH  Of  "reio^l  EprOUCNOY  = *,F14,7) 
IP  ( IiiOlbL  ,rO,  DGO  TO  20i 
PRINT  20«. 

202  tOPMAT(/,»  NO  AMPKNT  Nrisr  FIELI^  PRFSE’T,*) 
uO  lo  212 

2^1  HPINT  203 

r'>5  t''RMAT(/,»  AVRIEfjT  NOISE  FIELD  IS  CONSl  JERP"  , * ) 

?i2  ir( iHMr.ut.,EO,  1 )Go  TO  2r»- 
RWINT  204 

i-'.H  fORMAT(/,'  OMNIOlP'-CTIOr.Al  ELEMENTS  ARE  USE-*,*) 
oO  To  20b 


i'hiS  PAUi,  IS  BiiSI  QUALITY  PRACI1CABI4 
lX>rY  /UfciilsHl£D  TO  DQO  - 
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b7  ^0b  ?07 

bB  217  Format  (/»*  OIRECTIONAI.  ELf^ENTS  ARP  USE'^.*) 

b9  216  IFdNUrx.EO.DGO  TO  21P 

60  PRINT  POM 

61  20ft  F0PMAT</»»  GEMEPALIZEP  directivity  index  is  MAXlMi7fOd) 

o2  00  TO  ?09 

63  ilO  PRINT  PIJ  lALPMAOtfiETAUrGA'^'^AO 

64  211  FOPMAT</f»  DIRECTIVITY  INDEX  IS  MaXTWJZPD  WtTM  RESPeCT  TO  TTRECTIO 

66  IN  COSINES  (’lEia. 7, 

66  il9  1F(1ST8EP,E0,1»G0  TO  213 

67  PPirJT  214 

66  214  F1RMA7{/»»  the  array  IS  U*:STFrREp,n 

69  00  TO  21b 

70  213  PRINT  ?l6»ALPHAS,nrTAS»GAfc'<AS 

71  216  F0RMAT(/»*  the  array  IS  STEERED  WITH  ni^EfTtO*'  COSTt'ES  (*,El2.7, 

72  IS  • »tl2.7»  • » • »E12,7#  MM 

73  21b  Nr»0I.NTrl^*NDIV**242 

74  POINT  217»NP0INT,NP1v 

7b  217  fr*R^.AT(/,*  NU*mef/  IF  POINTS  CHOSEN  S'^ACF  z Mm»* 

7b  IS  CF  EACH  FACE  OF  ICOSAHEDRON  =MTS»MM 

77  iF(fJUlV.uE.l)GO  TO  3 

76  PRINT  2»nOIV 

79  2 FOPNaT(//»  ♦♦•••  ILLEGAi  PARAMETER  IN  S 'PRO  iTTME  RELAX,  NOTV  s M 

60  1 IS#*  IMPEl^lATE  RETURN  TO  ‘••AIN  "HOGbA",*) 

61,  Hf’jifHH 

62  3 lEtri.Gt  ,NRELAX)G0  TO  15 

63  PRINT  anal #n,nreuay 

64  2021  FORMAT(///*  NRELAX  CANNQT  RE  LARGER  THA  ' m,  ••  z • # 1 10 , » * oEl AX  z • 

bb  l#lin,*  RtTijRN  TO  haIN  PPOORAM,*) 

6b  RETURN' 

67  Ij  iFdSTELR.NE.nOO  TO  16 

66  SIMTsALPHAS*ALPHAS^PET,\S*  »ETASt6A'‘*'*AS*Ga»''»»S 

69  lF{SINlT,GE,,q99qo,ANn,sr  IT,lE,1,0«''«1)GO  16 

90  POINT  1 7 

91  17  FOPMaT(////*  array  is  SlEF'»En  PUT  THE  S' *M  Or  equarFS  OF  r'TPECTION 

92  ILOSInES  is  not  0i,E,*»/.*  return  To  MAIN  PRO^^RAH,#) 

93  RETURN 

94  16  IFdfJOEX.NEtDGO  TO  1ft 

95  SlNlTSALPHAO*ALPHAn^RETrO»PETAOAGAMMAn6'iAMM#.fi 

9b  1F(SINIT,GE. ,99999, AN;i*sr'lT,LFil,n00nllGO  tO  16 

97  PRINT  19 

96  19  FOHMaT(////*  directivity  INDEX  1$  TO  RE  MAXtHIZED  OUT  THE  SlD'  OF  S 

99  IWUARES  OF  DIRECTION  CoSIRES  js  hqj  ONE.*#/##  RETURN  TO  main  PROSRA 

100  2M,*) 

101  RETURN 

102  18  IF(N,rif  .NRELAX)GO  TO  2P 

103  IFdCOEF.LE.OdCOEESl 

104  lFdCOFF,GT.N)  ICOEFSN 

10b  GO  TO  21 

106  20  ICOEFSI 

107  21  PI2LAM:2.*3#14159265  /WAVLEN 

106  VSAVEs-1# 

109  nrelxi=nrelax-i 

110  SINITSO. 

111  lEdSTe£H,NE.9>G0  TO  4 

112  ALPHASsO, 

113  oFTASSD# 
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114 

lib 

116 

117 

llti 

119 

lao 

U1 

1^2 

UJ 

U4 

i2i> 

12b 

127 

1^6 

129 

IJO 

131 

132 

133 

134 
I3b 

136 

137 
13H 
139 
mn 

141 

14? 

143 

144 

145 

146 

147 
14(* 

149 

IbO 

Ibl 

lb2 

153 

154 

155 

156 

157 
15H 
159 
IbO 
I6l 
lb? 
163 
lb4 

165 

Ibb 

lb7 

Ibd 

lo9 

170 


wAMMASrD, 

4 uO  1 I:l*ri 

5Init=«;imt^ars(real<coef<I)))  > arS(ai  *Ar(corp^(T))) 

I UL  ( 1 1 sv  ( I ) •alphas* Y ( T ) •betas*’  ( I ) *ftAM»<AS 
IP(SINI r,GE.,noi)GO  TO  6 
Uf'  8 Irltfr 

8 C6fF(nsCMPLX(l.»0,) 

b ir{L00MA,GE,l)60  TO  q 
LOOP^Xr  25*((M/NHFLXl)*l) 

9 iF(N,EO.t.»ELAX)LOOPMXsJ 

ir(  ,ro,r)60  TO  l? 

PRINT  1718 

1718  KOR.«<AT  ( iMl  r • INITIAL  EXCITATIONS  ,,,  FLEMF^t  MNMPF9  <^’TTH  PFaL  AND 
ll**AOl(4AHt  PAHTS  of  ExCITATIONSM 
PPINT  J717,fl,COFF(I),Irl,N; 

1717  FORVaK  4(  I4*  * = ( • • , » ,ri?,7»  M • * t X)  ) 

I?  uOIVlSf  JlV-1 
JlV-? 

C 

c •••••*•«••••«•••••••«••«••«•••••••«*•• 

C uFGI..  r.rtU'iP  COOROInATF  PELAXATION  LOOP 


C 

uO  2r>0''  LOOP=l  »LrOPMX 
IF { IPPlNT.EO.riGo  TO  ?7 
CALL  FIHST 
27  ITlNDXr? 

uO  5 Irl»NRELAX 
un  5 J=1 »NPFlax 

A«^AIN<  I » J)=CMPLX(r,  »0,  ) 

5 A<IDr.  1 1 » J)=CmPLX(#',  »o, ) 
wO  7 l=2»NRtLAX 
r>r'X(  I )r  (L'^OP-l  I^npElxI  *1-1 
lNt>X(  I J=IM0X(  n-r.*(  ( 1NDX(  I )-l  )/N) 

ITslNOXd) 

uFfDscOtFdl) 

7 LPtFdl  )sC,N^PLX(0,,n,) 

C 

C COLLECT  THE  VEPTICFS  {>F  ThF  ICOSAmFORON 

C 

call  ArHAYS<NRELAX,AMAIr,,ASIDF»X,Y,7fnL»nTRrCT,COFF.P»l »l?»riny,M) 
IF  (NijIV'.tG,  1 )GO  TO  lO^'C 
C 

c collect  the  edges  of  the  TCOSAMEDPOM 

c 

CO  2no  J=I .NDIVI 
(‘rf.OlV-J 

^d  * 1 );j*p(  1 ( 1)  *K«Pd  *3) 

C(?»  l)rj*p<2»  d*K*P(?,3) 

3d)zJ*P(  3,  n *K*P(  3»3)  * 

Cd  »?)=J«Pd  (2)  *K*P(1 ,3) 

C(2i2):j«P<2f 2)*K*p(2f 3) 
u(3»2)=J*P<3»2)*k*P( 3»3) 

Ud  f3):J*P(l*4)*K*P<l»3) 

««(?»3)rJ*P(2»41*K*P(?»3) 

U(3»3)=J*P(3»4)*k*P( 3,3) 

G/d4)=J*Pd,4)*K*Pd»?) 


f 


) 
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171 

172 

173 

174 

175 
l7o 

177 

178 

179 
IriO 
101 
102 
103 
lft4 

105 

106 
107 
loO 
109 
19ft 

191 

19? 

193 

194 

195 
I9b 

197 

198 

199 
£0r 
201 
202 
C03 

204 

205 
^Ob 

207 

208 

209 

210 
211 
212 

213 

214 

^15 

2lb 

217 

210 

219 

220 
221 
222 

223 

224 

225 
22b 
227 


w(2«4)=J«P(2»4)>k«P(2»2) 

b(3»4)zJ*P(3»4)^K«P(3*2) 

Q < 1 » 5 > s J«P ( 1 • 4 ) ♦K^P ( 1 f 5) 

U(2*5):J4P(2*4)4-K*P(2»5) 

9(3»5):J4P<3*4)4K«P(3«5) 
w(l •b)sJ*P(l»4)>K*P(l»8) 

9<?>blsJ*P(2»4)^K*P(2*b) 

(i(3*6)=J*P(3»4)^K«PI3«6) 

call  APKAY5(NRELAX» Amain# ASI0r»XfV,Z*0L«nTPPCT»C0rPi0»5,6,lN0X*N) 
200  continue 

1F<ND1V.EQ,2)60  TO  lOftC 
C 

C COLLECT  TMF  FACES  CF  THt  tCOSAHEOPON 
C 

UO  3«''  J=1»N0IV2 

uO  300  Ksl,'J0lV2 

IFU  J4IC  ) .GT.NOIVl  )f;0  TO  300 

».=r0IV-J»K 

U(  1 * DzJftPU*  1)>K*P(  1 *2}  ♦L •PC  1*3) 

(»<2*  nrj«P(2»n^K«n<2>2)^L*P(2*3) 
s.(  3.l  )=j*P(3*l  )4K*P(3*2)>L*r(3#3) 

' «•(  1 i2)rj«n(  l*?)>K«r(  1 «3)4l  *0(1*4) 

u(2*2):J*P(2*?)'^K«ni2,3)4L«P(?*4) 

«(3*2)=J*P{3»2)*4*rC3,3)*L*r(3*9) 
o(  1 *3)rJ*P(  I »?)*K*P(  1 »5)*L*P(  1 ***) 
wl2*3)rJ*P(2*?)*K*Pi2*5)*L*P(?***) 
(.(3*3»rJ*P(3*2)*K*n(3,5)*l.*P(3*‘*) 
u( 1 *4)7J*P( t »4)*k«P<  t iS)*L«P( 1 *0) 

^■(2*4}:J«P(2«4)*k*p(2*5)4>L  •P(?*6) 
^(3*4)=J*P(3*4)*K*P(3,5)*L«P(3*61 

call  AFHAySCNRELAX,AMAlN»A5I0r#»*Y*?*'^L*nTP»C^*COfF,3#S,4,lM0X#N) 
30ft  C0^^T^4l.t 

c 

c polish  the  'MATRICES 

C 

1010  iTINiJAsl 

IFilMOf  X.UE.DGO  TO  no? 

^(1  * 1 )sAlPhAO 
y(2«  DztltTAD 
u(3*l):GAMMAU 

call  ARRAYS(NPELAX,AMAlr.#ASIDP*X*Y#Z*nL»DTPrC'^*COEF*0#l  #1  ,T:ir'X*N) 
iro?  LO  IftOl  1=1*NRELAX 

Jur.':ReAL{AMAIN(l#l) ) 

4‘'AlU(I#I)=CMPLX(mjM#ft,^) 

aiP:HEAL(ASIOE(l*n> 

AS  IDE  ( 1 * 1 ) ZCMPLX  ( 01  IN,  ft  , r ) 

LO  1001  JsI#NPELAX 
ASIJE(I*J>=AStOe(I,J)*APArN(I,J) 

AMAirK  J,1):C0NJG(AVA1N(I* J) ) 
irOl  ASIDEC J*I)=COMJG(ASlOE(Ir J) ) 

C 

c DO  the  eigenproblem  fop  the  group  of  COORPT^iATES 
C 

CALL  F rNISH(COEF#RATlO#NPrLAX, AMAIN# ASlUE#V,VX*0*WK*INnX#N) 
lF(WK(l),6T,0,,ANn,VSAVE,LE,WK(l))S0  TO  lR9a 
uO  1990  Iz2*NRELAX 
llrlNUXd) 
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i2a  cofp(in=oc(i) 

229  toO  TO  ?00l 

^30  19‘}9  lP(IPHINT,EO,0»6o  TO  19^7 

^31  CALL  SrC0»4D<HTIMr ) 

232  HPINT  1719,WK<1).L00P,HTI'^C 

^33  1719  f-OHMAT(//»  NEW  EXftTATlONS  ,,,LARAr<;T  ElGrN-ALUE  s •»F14,9»« 

234  1» •NFLAKATION  ITERATIOM  NUMBER  = ST^**  ...FiAPSEr  TIME  r ••2a6» 

«j3b  ^.'SECM 

23b  uO  172<*  Is2»NRELAX 

237  ITslNOXd) 

23b  1720  UIRECTd  isCOEFdl) 

239  IFdPUNCH.NE.DOO  TO  17?2 

240  KUf.CM  l723*(IN0Xd),|)IRtCT(t),l5?,NarLA<) 

241  17.?3  »«04-MAT(3dU,2Fll,6)  ) 

242  1722  PRINT  1 71 7,  dNOX  (1 ) * OIRK 7 d ) , 1=2 ,MRFL A»  ) 

243  1417  IF  (A»*Sd.-VSAVE/*K  d ) > .1  e .ARS(rPSI  dOO  fO  2-1'' 

244  2070  \<<;AVt=VAtU 

24b  C 

246  C •«••*«••••••••••••••«••••«••••**•••• 

247  c cNO  CiHOOP  coordinate  Rf|  AXATJON  LOOP 

248  C 

<49  c 

2b0  10  (LOOP'^A.ro,  I >00  TO  2010 

2bl  PRINT  2003 

2b2  2r''3  bOPMATC  MAXr^UM  fr'MuFR  QO  RFLAXATTONS  oovp-  OTFD  hTTHPUT  COrn/EPGFM 

2b3  icr**) 

2b4  oO  TO  POlO 

2Db  20^1  ppjiJT  2002 

2bb  2r''2  FOWMaT(*  ThE  OROuP  COf^RcTATE  bELATATTO’  mF^HOO  MAS  STOPPor  YIFlDI 

<b7  lf,r,  STEAOILT  I'^PHCVrO  RFSUt.TS,  CONyFROFNCF  TP  ASSUPFr.M 

2bb  2010  lOOPMX=lOOP 

2b9  lF(N,rit,r.RrLAXJGO  TO  157 

toC  uO  1 70  Kst  *N 

2bl  1 70  uC(K)rCJLF(K) 

262  1S7  ITCOfFrl 

<o3  ITJCslrOAd) 

264  ‘j7  hATrl0,*«L06lO(RATr0#VdTJC> ) 

2ob  PRINT  IOO.RAT»LOOP'^X,V<  ITJC) 

206  ion  FORMATdHl#/»*  GE.f4rRALl?£‘“  01  (DR)  * d Py , • - • *F2n,9 » / • • PU  *RFB  OF 

267  1 relaxation  iterations  r *»Tt7,/,*  TR"f  F tG^NVAL' 'E • » 1 7V , * r • » 

2b8  2820, <j/) 

269  b=0, 

2 70  LO  45  -=l*N  4* 

271  4b  bsS4REAUtC0eF(K)  )**2  + AIvA0(C0rr(|(>),*2 

272  lF(S.0fT,l.E-2'')G0  TO  44 

273  10  40  lsl,N 

<74  *0  UlReCTd)=CMPLX(o,,0,) 

27b  bO  TO  103 

276  44  brS9RT(i) 

277  UO  50  KSl»f4 

276  . COEF(K>=C0FF(K)/s 

279  AMAGsCAdS(COEF(K) ) 

260  lF(AMAG.LE,l,E-37)r,0  TO  40 

261  APHA=ATAN2(AIMA6<C0EF(K ) ) »REAL (COfF (« ) ) ) 

262  bO  TO  50 

263  49  APHAsO, 

264  50  UTRCCT(K):CNPLX(ama6,APhA1 
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206 

<67 

208 

209 

291 

292 

293 
299 

295 

296 

297 
296 

299 

300 

301 

302 

303 

p04 

305 

306 

307 
306 

309 

310 

311 

312 


^PWT»S 


103  PRINT  102 

102  KORMAT(»  OPTIMUM  CtEMCNT  EXCITATIONS  ARE#,/,*  rtT.  N0.*,t6X, 
1 • magnitude • * 23X  » 'PHASE • ) 

PRINT  l01»(K,niRECT(K)»Ksl*N) 

ni  KOPMaT(11O»2F30,9) 

1F<ITC0EF.E0,IC0EF)G0  to  m 
ITjCsIfiDxdl-ITCOEF 
iF(ITUC«LE,0)60  TO  56 
ITCOtFsITCOEF^l 

t/MAXsAhS<REAL(VX(1,lTJC)))4A8S(AlvA6(VX(l,Itjr))) 
lF(VHAX,GT,l,E-25)r,0  TO  l*>0 
uO  131  J:2»NHELAx 
oJslNDX(J) 

131  COEF(JJl=VX(J»lTjC) 

COEF(ri)=vxn#lTJC) 

00  TO  S7 

150  00  160  j=2,NRrLAK 
JJslNOX ( J) 

loO  COEF(JJ)sVX(J*ITjf )/vX(t • JTJC) 

COEF(N)=CMPLX(l,,r, ) 
toO  To 

56  ITJCsPOXdJ 

hATlOS^’ATIOAVUTJC) 

IFdCOf-F.E'i.DGO  Tft  SO 
DO  60  t^l*N 
oO  COEFCK)SDC<K) 

SP  KFTUHN 
tNO 
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STKtIT*HLLAX(l) .AHkAYS 

Sl»P»0UTlf4F  ARRAYS  wRTTTtN  RY  ROY  L.  ST'^F.TT,  ►'AVAL  UNO^RWATER  SYSTEMS 
CFNTLRi  new  LONDON  LAH0RAT0RY»  NFw  LONDON,  10  AU6UST  1977 


1 

2 

3 

4 
b 
b 
7 
b 
9 

1ft 

11 

12 

li 

14 

lb 

lb 

17 

lb 

19 

2ft 

dl 

22 

23 

24 

2b 

2b 

27 

2b 

29 

3ft 

31 

32 

33 

34 
3b 
3b 
37 
3b 
39 
4ft 

41 

42 

43 

44 
4b 
4h 
47 
4b 
49 
bft 
bl 
b2 
b3 
b4 
bb 
bb 


SUHHO'  TIME  ARRAYS  FILLS  THF  TWO  MATRlCtS  Far  THF  p'IQFmproHlE** 

ffRELAx  z 1 THE.  NUMHEP  rF  ELEMENT  CnR‘VF*iT^  TO  ‘RFLAY* 

AvAIN  z MATRIX  OF  QUADRATIC  FORM  FOR  MAiml-^^F  RFOTON 
ASIDE  z matrix  of  QUADRATIf  FORM  FOP  SIOFL^RP  RFOION 

x<i)  = X coordinate  of  I-TH  elemfnt 

Y(l)  s Y coordinate  of  I-TH  FLEmfnT 

7(1)  s 2 coordinate  of  I-TM  ELEm^nT 

'^L  = rOKK  AKF.A,  MFRE 
nio^CT  = WORK  Area#  HEPE 

COEF  z MRFStMT  VEfTOH  OF  El  E^FNT  CUPRF*'TS 
0 = APkAY  of  POirjTS  TO  ®E  POODfcTtr  O'lTO  ^HF  SPhFRE 
L1»L2  z loop  limits  DEPt'-’niNO  UPON  WUE'^F  T «F  POINTS  O I IF 
Tf.OX  = UiniCFS  OF  ELFMtNTS  RELAXFO  IN  ThIS  sTFP 
N s total  number  oF  rLPve-NTS  IN  TM**  AR^AY 

bUi’ROUTlf.t  arrays  (»  RELAX  , A ma  IN,  ASIDE » Y • Y , 7,  rL#  OIPpCT  ,COFF , 0,  Ll  »L?» 
1 TN()X,f.) 

ul'- EuSl  ON  <(n,Y(l),2(1  ) ,''L(  1>  lOt  **  I?)  fOoSOrTfS)  ,ST^’ROT(S) , IMrX(  i ) 
CO-^PL^  > MMAlN(NMtLAX,l),A«inE(NRF|  AX,!)  » p lOrCT  ( 1 ) , WTC  # COFF  ( I ) 

COf  MON  ir*OlSE»  I PHC ‘.E • ^AVL^N,  INDEX,  alpha  ', RF»  Aft  ,GAM‘’A0»  T COFF  , 

I IbTFEH,  ALPHAS,  nn  AS,  OA'^MAS,  NOT  V»'.CftP‘X,^PSI,  IPBINT,  IpuNCM 

Cft*  ><oi</ulR/lDlRl,  i-'ln^ 

CO*  M0N/£XTRA/PI2LA'-»  ITIfiDY 
uOUBLt  KHECISION  Pl'M 

DATA  CiSHOT/l.,.3ftftO I h9S4,-,Rft90U9‘^4,-,Ar*o-iSQ94,,3ft9ft  1^00(4/ 

LATA  SII*ROT/0, , .QStftSbSlb,  ,Sfl778S?S?,-,bP77**S^S2,-.9bl 
uATA  IiriKl,lOtH2/ft,0/ 

LATA  CrSl?/-l.E-6/ 

DO  3ftft''  L = 1»L2 

A = 1 ./SORT (Q(1,L) *0(1 ,L)40(2,L)*0(R»U)  + Q(3,L^  *^<3,1) ) 

alphis^ ( 1 ,L) WW 
prTAlsr.(^,L)*W 
0Amma:<M3,L)*W 

hOTAIE  lower  half  sphere  for  sYMmfTOY  R^ASO^'S 
IPlGAMs  A,ST,tPSI?)r,0  TO  3'' 

•/Ts-AuPHl*C0SR0T(4)*nETAl*SlNR0Tf4) 

orTAlS.rtLTAl*C0SR0T(4>-ALPHl*SlNRrT(4) 

ALPHlr*! 

j 

^0  jO  3ftOft  M=l »Ll 

AsALPHl*COSROT(M)*qETAUSriPOTCM) 
b=HFTAl#COSROT(M) ♦aLPM1*SINP0T(M) 

CfOAmma 

ALPHASA 

oFTA  SP 

lOIRslTlNDX 

lF(lNOrX,eQ,l)GO  TO  4 

CALL  REG10N(A,B,C»IDTR) 

1F(1UIR.GE,1.AND,I0IR.LE.')G0  TO  « 

PRINT  1»IDIR»ALPhA, BETA, gamma 
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b7 

b6 

59 

00 
bl 
t>2 
o3 

d4 

bb 

Ob 

o7 

bH 

o9 

70 

71 

72 

n 

74 

7b 

fb 

77 

7b 

79 

bO 

01 
h2 
o3 

64 

6b 

6b 

67 

Ub 

69 

9^' 

91 

92 

93 

94 
9b 
9b 
97 
96 
99 

ion 

101 

102 

103 

104 
lOb 
106 
107 
lUA 

109 

110 
111 
112 
113 


1 FO«MAT(//»  ••••••  •iUB'^OliTtNE  Pr^lON  GIV'-S  •••  TULrOAL  PPTUPN.V 

I • IniR  s rop  DlRCfTION  C0S!‘'F:S  alpha  5 •» 

* n5,A»*  ftFT^  s ••CIS.P**  GAMMA  s •,frt9.6) 

bO  TO  3000 

b lP(ll)lK.ta,3)GO  TO  30oe 
4 lF(lulP.LQ.2)1DlR?rl0lP2^1 
ir(nZP«L0.1U0lPt:2f)ZPl<»t 
-Tsl, 

IFUNOlSt,EO,")GO  TO  10 

AsAUPHA 

orPETA 

CsGAmMa 

ikT=0i<i0tbL(A*H«C) 
lF(><T.r.t.o,)bO  To  9 
POINT  7#<tT. aloha, BFTA.GAMVA 

7 ►04‘'<AT(//»  ••••••  SUHMOoTTNE  DNOX^E  GIVHS  a*  TLLrGAL  “rTUPr,*/ 

1 * fNuISt  2 *,F15,8*»  FOP  OlPFrTION  COFT* ALPHA  r «, 

C EIS.A**  BfTa  r •»E15,G»»  GAVMA  s •,Ft5,G) 

bO  To  'OOn 
9 i>JsSoUT{mT) 

10  1F< IPH04L,E0,O)6c  TO  ?0 
A:ALPHA 

u=PETA 

CsGA^ma 

isr. 

call  PmONE(A,P,C,OIPECT, 
bo  To  ?i 

20  uO  ll  U1,N 

11  DIf»ECT<I)=CMPUX(l,**'',0) 

21  •.TCSC'«PLX(0,,rt,) 

on  2b  IS1»M 

uOf-s(K(l)*ALPMATY(T)*HfTA*7(l)*6AMMA-nL(I)  )*pT?LAM 
A=COS(PbH»*«T 

D2SIN(«  oV)«MT 

bI«ECT(l)sCMPLX(A,fl)*OlHECT{ X) 

2S  »,TC2k»TCTC0EFn)4Dl9ECTIl) 

«TC2s«rAL<bTCJ*PFAL<MTC)  ♦ A I '♦AG  (nTC ) ♦A  ImaG  , ^-rc) 

• TCsCOfrJblwTC) 

1F(IJU.E9,2)G0  TO  2O0P 
AMAlN(l*l):AHAXN<l,l)4WTC;* 
oO  100-  ls2,NPELAX 
lUlNDXd) 

MMA I N ( 1 , n CAJ4A  IN ( 1 , 1 ) *01  RFC T < T I ) bWTC 
oO  100''  j:I»NRCLAX 
JJsInOx ( J) 

irno  AMAlr4<l,J)=AMAIN(I,j)4C0NJG(DIRECTfII))*0TPrC^(JJ) 
bO  TO  3000 

2000  A9X0En  •l):ASIDE(l,l)<fWTC2 
UO  2001  l:2»NnELAX 
nslMOXd) 

ASmEd  dlsASlOEd,!  INDIRECT  (IT  UwTC 

uO  2001  jsX»NR£UAX 

JJslNDXtJ) 

2001  ASIOEd*  J):AS10E(X,J)4C0NjG<0IRECTdn  )«OTRrCT(JJ) 

3000  CONTlNUfc. 

RETURN 
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I T#mLLm»  ( n 


I c 

? c 

J c 

“ c 

t 

b C 

7 I 

h C 

C 

10  c 

n c 

u c 

13  c 

w c 

c 

lt> 

17  c 


1« 

1*^ 

^1 


i»"(«0uTlNf:  P’INISh  wRITTfN  »Oy  *'AVAL  UNHEPwAT^P  SYSTEMS 

CFnTEM,  LO'JUON  LA^OwATOQY,  NFw  L0**D''N,  rT  in  AUfiU<;T  1Q77 

<;uftHOi.TlNE  finish  solves  TmE  EIGFMPRORLEm  ^’jn  RETiJPN<i  TMF  COFFpirlFNTS 

cotF<i)  s THE  Excitation  coecfictfnt  of  the  t-tm  FLE»*F'fiT 
PATIO  = estimate  nF  THfc  PATIO  OF  POWER  IH  ▼h*'  MA1NL0*^F  REGION  TO 
POWER  If]  THE  SlO^’LOhr  REGION,  MTT^  TS  MOT  IN  OP, 

NmELAx  = NDMREH  OF  RELAXED  FlEMFmTS  ♦ ’ 

AVAIN  r MATRIX  OF  QilAnPATlr  EORm  FOP  MAikiLARe*  REGION 

AsioE  = matrix  oe  O'jAnHATrr  form  fop  sidfl-*rf  region 
UjOx  = INOICFS  OF  elements  RflAxeO  T»'  TmTS  STEP 
N s NhmIjER  of  ElE-’E'ITS  r THF  ARpAV 

ThF  RiMAINOER  OF  THE  AF6"mfntS  ARF  AS  ^FOMtRfo  SI 'PPOUT  J NT  PEMCL  H 

S'»|<ROiiTlNfc'  EINISM(COFF,PATIO,'''RFLAX*AMAlM,Aet^F»V,VX,n,WW,lNnx,N) 
uTMtfiSl  Jli  V(  1 ) ,U(  1 ) » VK  ( 1 ) » INOX  ( 1 ) 

c0‘  PLt  y AMAlN(NRrL^x,  1 ) , Astf^E  (MPpi  AX, ! ) ♦ vx  ( ‘ RELA X , 1 ) ,rorr{  1 ) 

Cn^  MoN/JlP/lUlWl  , I-'IU? 

' rzNRLLAX 

,opt=i 

lALL  Pf '^CLH<  AMAlfj,  aSI'F  »m,v,  VX,»*,MonT,c  * K,»innn) 

) 

IF  ( JJ.F  v<,n)G0  70  l^On 

1F( jj.i  j.nrFlaxigo  to  1*^ 


c7 


Ojjs  JRV  lax-jj 


il 

s? 

i3 

Jf> 

37 

3n 

3N 

^C 

m 

4^ 

43 

44 

45 
4t> 
47 
4ri 

49 

50 
bl 

53 

54 

ss 

Sb 


JOINT  rOOl.JjsJ 

{-oimati///’  xxxxy  n.itficai  ^'IEf ioulti'-s  i»  ''he  FioeMAHALYsis  prf 
ivP  jT  TJ'L  CO‘U»LETf  SOU>»TJON  OF  TOlip  PpOSLEv.,/ 

C • a/xXa  The  Hj-ST  of  THF  CO*'PUTrO  C''ETFTrTE»'TS  is  RETI'^»'EE,  •/ 

3»  xy<X>  TMfRE  Al,r  AT  SIS,*  LT^;EA  -Xv  •N''EPrN'TF.NT  RFTTFR  SfTS 

4 OF  COf  FUClENTS.*///) 

IS  1T=1 

V*‘AX=-1  ,t  ♦2'' 

?('  jrl,JJ 

VFlN13SuT,N1'AX)r,0  10 
w‘’AX  = Vf  J) 

ITrJ 

:n  tONTlNi't. 

INI  xai  = n 

n*'  ( l)=‘v  ^AX 

.<A7  lor  FI  OAT  < T0I(»?  > /f^OAT  ( IniRI  ) 

w''AXrA'*S(RFAL(  VX  (SIT)) ) ♦ a PS  ( A I^Ar,  f V)i  f I » IT)  S 
IF( ^A>  ,yT,  I ,F-^S)'‘.0  TO  S'' 

bO  50  J=l,N 

30  oOrnJ)sCOFF(J)*vXtl«IT) 

JO  31  jr^, NRELAX 
jjrlN[jy(d) 

31  COF F( JJ)rVX ( J, n ) 
oO  TO  1002 

‘,0  dO  bO  jri',NRL'LAX 
JJrPJDy(J) 

f>0  C0EF(JJ)=VX(J,IT)/VX(SIT) 
oO  TO  JOO? 
ino  HRIfIT  1001 


. J 

\ I 


n-  in 


n«  IS  B»SX  OUALIW  PBACSWABW 

n»«P^®?iil,c«nXODDO 


i 
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1 


b7 

b« 

b9 

t»Q 

bl 

62 

63 


lOOl  f^RMATi///*  H»XXX  NuvrMfAL  niPFIClH.TIES  T^*  TMC  T IfiENArMl  YSIS  PRf 
UP^r  7i<e  SOLUTION  nr  rOl  R PROPLE'^,  sorry,  CmARLCY,  •///) 

ir''2  ioi«is^. 

iniH2='' 

NFTURN 

£Nr 


IK  r)8«9 


SThL I T •KtLA* ( I 


1 c 

^ c 

L 

**  c 

i>  C 

^ C 

7 L 

« C 

^ C 

I*'  L 

n c 

13  C 

1**  t 

1*3  L 

L 

17  C 

IM  (. 

1^*  c 

c 

^’1  c 

L 

CJ  C 

4i'*  C 

c 

c 

t 

^0  L 

C 

c 

31  L 

32  (. 

33  c 

3**  C 

3t>  C 

3o  f 

37  C 

3m  C 

3‘>  C 

“♦I  C 

•*2 

•43  C 

•4«* 

•4b 

Hb 

•♦7  C 


•4b 

•49 

bO 

bl 

b2 

53 

b4 

bb 

bb 


b'iHROU7lf4F  PEtJCLh  WRITTEN  RY  ROY  L,  RT*<ETT,  NftVAL  UNOFPwATf.R  Sy5TF.ms 
cr^TCW*  NFW  LONDON  LAROPATOPY,  Nflk-  LONO'^N,  fj  17  PAY  IQTF 

'•.liHROrTlNF  PFNClh  computes  ALL  THE  Fl6fNVA|ijrS  ANO  ALL  TMp 
PIOf-lVtCTORS  OF  T''E  OEnEPAlIZFD  riRENP'^OBLrM 
A • 7 : LAMROA  • R « 7 

w»  EHt  A AriD  0 APF  HFRMITIAM  vaTRICFS  a iu  R IS  POSITIVE  OFFinITF, 

TtaS  ROUTINE  USES  TwO  PROGRAMS  OF  I^SL  INTERNATIONAL  ir 

A^'0  statistical  libraries, IMC,»  hOUSTO'N  TrXAS)  FOR  ThF  SOijjTTON  OF 
TMf  Or.ulNARY  HFR‘*ITIA^  FIOENPRORLF”  (^HAT  IS»  WHFN  TtlF  MFiJMlTlAN 
“ATiUv.  ri  IS  JUST  THE  iDf^TITY  PATPl*)f  ANY  ^FT  OF  PO'‘TI»;FS  FPP  SOLVl^’O 
TmIS  PHUMLEF  CArj  RE  SUhSTITDTEO  fF  CARE  TS  '|FFD, 

T.-I'a  POl;TINL  EXTfOS  TC  THE  ‘tFHMjTIAf  CaSF  tmF  ALOORTTMV  op 

»J,S,  vAh*TlN  ANO  J,H,  tPFOtlCTlON  OF  THF  SYM^TTPIc  FTGEN- 

PUO'^LfM  \,x  = LAVRDA  n*X  ANO  RFLATLD  P-^OMLEmS  TO  STAMp^pO  FORM*, 

•ii  •:l‘>,  'ATM.,  \/oi  , ii»  p,og-iin,  iof,«,  wnT-'4  TS  DFSiorrn  f-ip  the  rfaj 
CASE. 


CAoTIPf.  •««  ALL  Of  A A j ' Tm'’  ‘TRICT  LOw^P  TOI'N-^'LE  PF  !!  APf.'  ''rsTRpYFn, 

A,;4  z nfPMlTlAfj  •AT.'ICf'-  OF  TMF  F lOC-^WRrn,  COMFLFX  '-'ATiUOfS 

^ I <*Ff»SIOf  F'*  IN  X IM, 

Mr  r.  L OWOEP  OF  TH^  FIG^NPRORLFm 

TL  s TJL  ACT"AL  Pn^  nlf/E"SIOM  OF  ROTM  a am-  n, 

VALJt  = WILL  COt.TAl-j  The  EIGFMVAiUES,  <EAL  ARRAY  ''IMr*.|*^Ipv*Fn  A7  iEaST  *1, 
VfCTO*-  s will  COriTATM  THC  F IfiFNvrCTORS.  CO;  m“N  J OF  VFCTOR  COMTATMS 
THE  EIgFmVFCToR  COHmESPoNOIMG  TO  EIGENVALUF  VAL‘'f<j), 

C0»APLLX  AR^AY  DIMENSIOmeO  IV  * M, 

IV  S RJW  OIMFNSIO'I  OF  VECTOR,  **UST  HAV‘  TV.0F,N 
NO»M  = OPTION  ImC.TCaTOR, 

z o , COMPUTE  F1GLN\/A|UE5  ONLY, 
r I » COMPUTF  riGFNVALllFS  A**''  FIG*- NVFwTORS, 
ol  s mEal  array  pt«f:‘isio*  fp  aT  lfast  n.  wo-k  arfa, 
w*-  r ( LAL  array  nT'<r'iSio*"rn  at  least  3*n.  opk  area  'vhosf  co^-temts 

AHt  nESTROYFH,  ON  Nrp‘.«AL  RETMRfI,  CONTAINS  THE  NU^HFO  Op 

CJWRFCTLY  C0‘'PMTFL'  EIGENVECTORS,  ^ PPOf»  »FT»IR')#  r -1, 

% z h ilATEMFNT  ri'jMMFR  P-  THE  CALLING  RR^'G^'A^,  RETURN  TO  STATEMENT  S 

O'LT  IF  t4  IS  MOT  POSITIVE  definite.  rC'M  GE  CAUSED  MY  MOHMPP'C,  poROR) 

b'JhWoUT  ir.f  PEMCLM(  A»H,N,  I*  , value  » VECTOR*  I V » • OPT#  PL  » V*K  » $ ) 

uI'^Et^SlU..  A ( IN#  1 ) ,t>(  1*1,  I ) , vALUFM  ) , VFOT'^R  ( !»  , 1)  » OL  ( 1 ) » VK  { 1 ) 

LO'PlEX  A»H»vrCTOP,X 
uOUMLt  PKFCISION  2 

bFGiM  **AHTIN  ANL  hTLKINSO'  ROUTINF  RFOU':! 

ifin,gt,ojgo  to  S 

rjr»N 
GO  To  11 
b *>K(2)Z-1, 
uO  1 1 = 1 .rj 
LO  1 J=l»N 
x=C0rjJG(b(  I » j) ) 
lF(I,LO.l)GO  TO  4 
lPsI-1 


H-  I J 


I 

TR  5889 


M 


t>7 

uO  3 KrlP, 1 »-! 

bH 

3 

XzX»CONJb<H(I*K) )*n(J»K) 

b9 

4 

iri  TO  2 

toO 

inweAL  (x),le*o.)rfTii«in  n 

oi 

IF  ( AHS{Al*1AG<X)/f«f  AL(X)  » ,GC,l  ,F-4)RF:TtfH  . 1 1 

6? 

2zREAL(X}««2 

bJ 

4:2tAlVAO(X>*42 

o4 

uL(  ns9UKT(S0PT(Z)  ) 

bb 

60  To  1 

bb 

2 

0(J»n=A/DL(l) 

b7 

1 

tO'-TI*-*!  It 

oB 

ll 

UO  2''  1=1. fi 

uo  2''  J=I»M 

XzAd.Jl 

71 

K { I ,t( . noo  TO  14 

72 

10=1-1 

cO  it,  czio.l.-l 

74 

Lb 

A=>-M( I .A )«A< J.K ) 

7b 

14 

«<j.nrA/.bL<n 

7o 

COf.TIMifc 

77 

wO  IP  J=l  »fl 

n 

uO  3 ^ J =0 . ' i 

/S» 

a=A( l.j) 

nP 

K c I ,cf  . J)oO  TO  44 

(11 

i'*=l-l 

b2 

^0  3b  ► =1P. J.-l 

b3 

lb 

a=X-A(i-  .J)*C0MJb(P(I,>^)  ) 

04 

•*4 

lP(J.t'*.noO  TO  4^ 

bb 

IP=J-1 

Ob 

uO  3b  K=1P,1.-1 

07 

lb 

x=x-C0r  Jo( A( J.K ) •n( I ) ) 

OM 

*b 

A(  I.j)  = x/0L(I) 

o9 

IP 

CONTllii/c. 

90 

C 

Finn  T(tt  EiGrbyALUf-s/ricF:*  VFCTOPS 

91 

M 

C 

c 

92 

T=IEAL(A(I,I)) 

9J 

4n»n=c-^pLX(Y»o.) 

94 

uo  so  J=I ,0 

9b 

A(1>J)=A(J.I) 

96 

AiJ.Drtu  ur»(A(j.T} ) 

97 

(. 

I'^SL  Rrjri'.E 

90 

CALL  VCvTCH(A*U* IM*A) 

99 

|r<fioPT,LT,'',OR,NOPT,&T,?)NOPT=l 

100 

ir(M0PT.  -F,''.AND,TV.LT,fc)f  0PT=P 

101 

t 

i‘>SL  Rf  o!  ne 

102 

call  f.igch(a  »n,m''pt. V7li'F.vfctop, tv.wa, rr*. > 

103 

■».PsN 

104 

1F(  ItR.ot  .laHl^'asIrM-lPB 

lOb 

IF(M2.ot.l»G0  TO  Si 

106 

..K<2)=P. 

107 

RFTURN 

106 

•Jl 

AM  {2)=^‘2 

109 

lF<NOPT.tO,0,OR.»iOpT,FO,3)PFTijP^ 

110 

uO  9P  J=l,‘<2 

111 

2=0. 

112 

uo  flo  rsl.f/ 

113 

60 

2=/*PeAL(VrCT0Rl  1 , j)  )**i;AAtMAf;(VFrTOP(  I » J)  )••? 

B- 1 5 


i 


lU 


ll'*  Y=SQKr(^) 

Ub  wO  I = l»;4 

ilh  «1  vrCT0»<(l*J)=VfCTOMfI*J)/Y 

117  40  tO^T^N(^t 

Urt  C "mHTIN  and  ^.TLkIMSO*  hOtiTU.'F  ^FHaka 

11^#  DFi  jrl#v2 

12''  1 S'**  1 »-l 

1^1  x=vfcC7r*K<  I , j) 

122  IF(I.e  ;.f.)r>0  TO  71 

12J  10=1*1 

12**  jo  T''  * :iP»M 

I2b  70  AT*-CO/.s/onUK»I)  )#vtCTO>*(rf,j) 

UO  71  wrCTO«n»J)=X/DL(I) 

127  ,C  LOr.TlNlc. 

hFTkJW  J 
tf'? 


fcP-t  T»b 


XUIS  PAO£  IS  BEST  QUALITY  PRACXICABLI 
IBOd  OOPY  FWWISHfiD  TO  DDC  


J 

i 
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sthext*m: 

1 

i; 

i 

<*. 

5 

L 

7 

a 

9 

IS 

11 

12 

13 

19 

lb 

l(> 

17 

lit 

19 

^r 

21 

22 

29 

2b 

2b 

27 

2d 

29 


LAXdl 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


.RE.&ION 

bUPROUTlNF  region  written  BY  ROY  I..  STREIT,  NAVAL  UNDERWATER  SYSTfmS 
center#  NEW  LONUON  LABORATORY.  NEW  I ONDOn.  rT  10  AUGUST  1977 

suuhoutime  region  Divides  polar  sp/\ce  into  three  disjoint  regions 

AlPHA  = direction  COSINE  WITH  RESPECT  TO  THE  X-AXIS 

beta  = direction  COSINE  WITH  RESPECT  To  TmE  Y-AXIS 

GAHMA  = direction  COSINE  WITH  RESPECT  To  The  7-AXIS 

iriR  = t IF  DIRECTION  COSINES  Lir  IN  THE  WAINLOBE  REGION 

: ? IF  Direction  cosines  li^  in  the  sioelobe  region 

= A IF  Direction  cosines  lie  in  the  ignored  region 

SUHROUTINE  REGIOn(ALPHA.BETA. gamma. tDIR> 

•*****«••*•****•*•**••*•*•*•***•«*•*• 

OEGIN  code  for  specific  array  problem 

**t******m******* ********* *********** 

SET  COSINE  OF  82  DEGREES 
DATA  COSfa?/. 1391731009/ 

IDIR=2 

lF(AUS(bETA) .LT.COSe?) IUIR=l 

**•*«*• *•••********••«•* *••«•••*••* 

EtJO  code  for  SPECIFIC  ARRAY  PR0BlE“ 


return 

ENC 
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I 


I 

f 

I 

I 

I 


I 


SlHLlT*KtLAx<i) .DNOlbE 


1 

2 

i 

4 

b 

( b 

} 7 

I fc 

, 

, 

I 1 
12 
n 

14 
lb 
It) 
17 


C bHhROUTliJE  ONOISF  rtRiTTCN  F3Y  ROy  L,  mavAL  UNnrPWATFP  *;YSTFMS 

C CrMEM»  UE*)  LONOOH  LARORATORYi  Nfw  LONQ  )n,  ej  10  AU^^UST  1Q77 

C 

c riiNCTlON  ONOISE  C'^MPuTES  the  NOKF  level  I*  ANY  OIRECTIOm 

C 

c alpha  s direction  COfilNE  «^ITm  RESPECT  To  TwE  Y-AXTS 

C itfTA  = DIRECTION  CO<;lNE  -•ITM  RF«;PpCT  Tq  JuF.  Y-AXTS 

c f^AM-'.A  = DIRECTION  Cosine  /<ith  respect  To  Tj.e  z-axts 

c ONOISF  s NOISE  LEVFL  IN  T'^E  SPEClrlEP  OIRECtION 

c 

f I INC  T I On  'JMOlSE(ALPMA»BrTA,r,AHMA) 

C 

c 

c oEMN  COoE  FOF<  SPECIFIC  AORAY  PROPLE'-' 

c •••••«••«*•••••*«•••••••*««•«•*•••*«• 

c 

^j'lClSEzl, 


lo 

19 

in 

21 

22 

23 

24 


c *«••«•*••••*«•*»•«***••••  •«••**••• 

c t'lr  Curt  FOR  SPECIFIC  array  ppoblfm 

c «*•*«*•«******•**••*«*•«•••**•••••• 

c 

rETOKN 

L*n' 


.Pn»b  M.f'Mol.E 


I 


i 
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TR  5889 

I 

! 

STHtlT»MELA> ( 1 ) .PHONE 


1 C bUBROUTlHE  PHONE  WRITTEN  PY  ROY  L.  ST“EIT,  N»VAL  UNOERwATFR  SYSTEMS 

2 C center,  new  LONUON  LAHORATORT.  new  LONOHn,  fT  in  AUnilST  1977 

c 

9 C SUOROI'FINE  PHONE  COMPUTES  THE  RESPOF'SE  PATtF9N  OF  EACH  MYOnoPHONF. 

5 C the  USlH  Ml/ST  Take  care  not  to  change  The  y,  y,  and  7 arrays. 

• b c I 

7 c ALPHA  : niHECTION  COSINE  WITH  RESPECT  To  TuE  K-AXIS  ! 

8 C beta  r niRECTION  COSINE  WITH  RESPECT  To  TwE  Y-AXIS 

9 C gamma  : niRECTION  COSINE  WITH  RESPECT  To  ThE  7-AXIS 

15  C OIRECK  1 1 I . . . . DIRECTil.)  : THE  DIRECTIONAL  RESPONSE  OF  I:  HyoROpHONFS 

, 11  C IN  The  SPECIFIF-  OIRECTIO'I. 

12  c (note  that  The  response  m*y  pf  ro“PLEX) 

li  t N : TOIAL  NUMBER  OF  HYDROPHONES 

19  C X(I)  A COOROIIiATE  OF  1-TH  HYDROPHONE 

lb  C Y(I1  : Y coordinate  OF  I-TH  hydrophone 

lb  C /(ll  : A coordinate  of  I-Tw  HYDROPHONE 

17  L 

18  b'.IHHOUTlNE  PHONE(ALPha,hetA, GAMMA, DIRECr.H.v.Y.Z) 

19  C 

20  COAPlEY  OIRECTll)  ' 

21  uIMENSION  X(1),Y(1),Z(I)  I 

22  C i 

23  L •••••••••••••••••••••*»•«»•••»••*».««  1 

29  c uEGl!.  CODE  FOR  SHErlFIC  ARRAY  PRORLF*'  I 

2S  C 1 

2b  c : 

27  uO  10  I:1,N  I 

20  in  i,IHECI(I):CMPLX(i.n,n.o, 

29  c 

30  c ....... A.*.*. 

31  t L'll  CO!'u  FOR  SPECIFIC  ARRAY  PROBLF" 

32  C ................................... 

33  C 

39  HFTURN 

3b  c-ND 


. FIN 
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